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Q . Abstract 

CN ! We construct line shapes for the X(3872) that generahze the Flatte and zero-range hne shapes 

bJQI that have been considered previously. These line shapes are associated with scattering amplitudes 

that are exactly unitary for real values of the interaction parameters and can be derived from a 

renormalizable quantum field theory. The new line shapes can be used to discriminate between the 

alternative binding mechanisms in which the X(3872) is generated either dynamically by charm 

meson interactions or by a resonance near the D*^D^ threshold. If the resonance is identified 

' f— I 'l with the P-wave charmonium state Xci? the interaction parameters can be constrained by using 

Qh' charmonium phenomenology. We analyze data on the X(3872) and also data from the Belle and 

pi • Babar Collaborations on the invariant mass distribution of the charm mesons from the decay 

^ ! B ^ K + D*'^D^ up to 4000 MeV. Our analysis is compatible with the mechanism for the X(3872) 

being either a fine-tuning of charm meson interactions or the fine-tuning of the Xci to the D*^D^ 

threshold. In particular, the data do not exclude a separate x'ci resonance between the D*^D~ 

> ■ threshold and 4000 MeV. 

l> 

^ ■ PACS numbers: 12.38.-t, 12.39. St, 13.20.Gd, 14.40.Rt 
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I. INTRODUCTION 

The discovery of the X(3872) resonance by the Belle Collaboration in 2003 [1] marked 
the beginning of a new era in charmonium spectroscopy. About a dozen new cc mesons 
above the open charm threshold have been discovered and many of them have properties 
that seem incompatible with their identification as conventional charmonium states. This 
presents a serious challenge to our understanding of the cc sector of QCD. Of all the new 
cc mesons, the X(3872) is the one for which the most experimental information is avail- 
able. The preponderance of this information implies that the X(3872) is a charm meson 
molecule whose constituents are a superposition of D*^D^ and D^D*^. The structure of this 
molecule is remarkable, with the charm mesons almost always very well separated. However 
this identification of the X(3872) is not universally accepted within the field. The leading 
alternatives are the ^Pi charmonium state Xci (2-P) or a compact tetraquark cc meson. Even 
among those who lean toward its identification as a charm meson molecule, the remarkable 
structure of the X(3872) is not universally appreciated. 

The only experimental information that is necessary to make the identification of the 
X(3872) as a loosely-bound charm- meson molecule is the determination of its quantum 
numbers and the measurements of its mass. The quantum numbers of the X(3872) can be 
inferred to be 1"^"*" by combining the following information: 

• the observation of its decay into J/ip'y [2, 3] or iIj[2S)'j [4], which implies that it is 
even under charge conjugation, 

• analyses of the momentum distributions from its decay into J / ip n^ n^ , which imply 
that its spin and parity are 1+ or 2^ [5, 6], 

• either the observation of its decays into D^D^n^ [7], which disfavors 2^ because of 
angular-momentum suppression, or the observation of its decay into 'ip{2S) 7 [4], which 
disfavors 2~ because of multipole suppression. 

A recent analysis of decays into J/ip tt+tt^tt'' by the Babar Collaboration favors the quantum 
numbers 2 ^, but does not exclude 1++ [8]. In the absence of definitive evidence to the 
contrary, we will assume that the quantum numbers of the X(3872) are l^"*". 

The mass of the X(3872) can be determined most directly by measurements in the 
J/ip'K~^7i~ decay channel. As pointed out in Ref. [9], measurements in the D^D^ti^ decay 
channel are biased by the associated threshold enhancement just above the D*^D^ threshold. 
As pointed out in Ref. [10], measurements in the D*^D^ decay channel are biased by the 
analysis procedure in which D^ir^ with invariant mass near the mass of D*^ is constrained 
to have invariant mass exactly equal to M^,q. This procedure assigns an energy above the 
D*^D^ threshold to a D^D^tt^ event whose energy is below the threshold. Using the most 
recent measurements of the mass of the X(3872) in the J/ipir^ii' decay channel by the 
Belle, CDF, Babar, and DO Collaborations [11-14], the combined average for the position 
of the X(3872) resonance relative to the D*^D^ threshold is 

Mx - (M,o + Mo) = -0.42 ± 0.39 MeV, (1) 

where M*o and Mq are the masses of D*^ and D^. 

The reason the quantum numbers l^"'" and the measurement in Eq. (1) are sufficient 
to determine the nature of the X(3872) is that quantum mechanics implies that an S-wave 



resonance whose energy is sufficiently close to the threshold has universal properties that are 
determined by its energy. If the energy is below the threshold, the state is a bound molecule 
consisting of pairs of particles that are almost always very well separated. One of the 
universal predictions is a relation between the mean-square separation of the constituents 
and the binding energy Ex'- (^^)x = (4/i-Ex)~^j where n is the reduced mass. In the 
case of the X(3872), the quantum numbers 1"'""'" imply that there is an S-wave coupling 
to D*^D^. The tiny energy relative to the D*^D^ threshold implies that it is a resonant 
coupling. Thus the X(3872) is an S-wave threshold resonance. The binding energy given by 
Eq. (1) implies that the root-mean-square separation of the charm mesons in the X(3872) is 
{r'^)x = 4.9l}^3^ f'^- This huge separation of the charm mesons makes the X(3872) a truly 
remarkable hadron. 

One should distinguish between the nature of the X(3872), which is a weakly-bound 
charm meson molecule, and its origin, which has to do with the binding mechanism for the 
molecule. Identifying the origin of the X(3872) is crucial to understanding its implications 
for the other new cc mesons above the DD threshold. Universality is a double-edged sword. 
While it allows the nature of the X(3872) to be determined unambiguously from limited 
experimental information, universality makes it more difficult to identify the origin of the 
state. There are two primary candidates for the binding mechanism of the X(3872): 

• dynamical. The X(3872) could be generated dynamically by the interactions between 
the charm mesons. The interaction between D*^ and D" in the isospin-0 1"'"'" channel 
could be tuned to near the critical strength for the formation of a bound state. 

• resonance . The X(3872) could be generated by an isospin-0 resonance whose energy 
is tuned to near the D*^D^ threshold. An obvious candidate for the resonance is the 
Xci(2P), but it could also be any other type of cc meson with the appropriate quantum 
numbers, such as a compact tetraquark. 

One set of clues to the origin of the X(3872) is its decay pattern. There are 6 decay 
modes that have been observed thus far: J /ipn^n~ ^ J /ipTi^-K^-K^ [3], J /ip'~i [3, 4], ip{2S)'~^ [4] 
D^D^TT^ [7], and D^D^'j [15, 16]. In order to exploit this information, one would need to 
understand the pattern of branching fractions that follow from each of the binding mecha- 
nisms. There is a well-developed phenomenology for decays of charmonium states, but there 
is no corresponding phenomenology for decays of compact tetraquark cc mesons or for the 
inelastic scattering of charm meson pairs. This makes it difficult to constrain the origin of 
the X(3872) using the observed branching ratios. 

Another class of observables that can provide clues to the origin of the X(3872) is the 
energy distribution or line shape in specific decay channels. At energies much less than 
8 MeV from the D*^D^ threshold, the line shapes in these and other channels are universal 
in the sense that they are determined only by the binding energy and width of the X(3872) 
resonance [9]. The universal line shapes in the D^D^n^ and D^D^'-f channels, which receive 
contributions from the decay of a constituent -D*° or D*^, are different from those in other 
decay channels, such as J/ip tt+tt^. Braaten and Lu presented simple analytic expressions for 
the universal line shapes that take into account the width of the D*^ and inelastic scattering 
channels for D*^D^ [9]. The accuracy of these line shapes in the D*^D^ threshold region 
was verified in Ref. [17]. 

There have been several theoretical analyses of the line shapes for X(3872) produced by 
the decays B ^ K + X. The specific decay channels for which the line shapes have been 



measured by the Belle and Babar collaborations are J/ip-K^Ti' [11-14, 18, 19], D^D^n^ [7], 
and D*^D^ [15, 16]. Hanhart et al. used a Flatte parameterization of the line shapes to 
analyze the data on the J/ipTc'^'K^ and D^D^ir^ channels [20]. They concluded that the 
data favored the X(3872) being a virtual state with energy above the D*^D^ threshold. 
Their analysis was biased towards such a result, because they assumed that the line shape 
in D^D^Ti^ vanishes below the D*^D^ threshold, thus ignoring any contributions from a 
resonance below the D*^D^ threshold. Braaten and Lu used the universal line shapes to 
analyze data on the J/ipTi^T^" and D^D^n^ channels [9]. They concluded that the data 
favored the X(3872) being a bound state with energy below the D*^D^ threshold. Zhang, 
Meng, and Zheng [21] followed Ref. [20] in using the Flatte line shapes and ignoring D^D^n^ 
and D^D^'-f events from a resonance below the D*^D^ threshold. They concluded that 
X(3872) can be identified as a ^Pi charmonium state that is strongly distorted by coupled- 
channel effects. Braaten and Stapleton used the universal line shapes to analyze data on 
the J /ip 71^71^ and D*^D^ channels [22]. They pointed out that the analysis procedure 
for D*^D^ in Refs. [15, 16] assigns an energy above the D*^D^ threshold to a D^D^n^ or 
D^D^'j event from decay of a resonance just below the D*^D^ threshold. When this effect is 
taken into account, the analysis of the D*^D^ data favors the X(3872) being a bound state 
with energy below the D*^D^ threshold. Kalashnikova and Nefediev [23] followed Ref. [20] 
in using the Flatte line shapes and ignoring D'^D^tt^ and D^D^'y events from a resonance 
below the D*^D^ threshold. They concluded that the Babar data prefers the X(3872) to 
be a virtual state with a small ^Pi charmonium component while the Belle data prefers a 
bound state with a ^Pi charmonium component of about 30%. 

Since the lines shapes very close to the D*^D^ threshold are universal, information about 
the origin of the X(3872) can only come from the energy distributions away from the thresh- 
old. The line shapes outside the universal region depend on the binding mechanism. The 
Flatte line shapes used in Refs. [20, 21, 23] can be derived from the assumption that the 
charm mesons scatter only through their couplings to a resonance with isospin 0. These 
would be the appropriate line shapes if the X(3872) arises from a tuning of the energy of a 
resonance. Alternative line shapes have been derived by Braaten and Lu under the assump- 
tion that the charm mesons scatter through zero-range interactions between the coupled 
channels consisting of neutral charm mesons and charged charm mesons [10]. These would 
be the appropriate line shapes if the X(3872) arises dynamically from interactions between 
the charm mesons. The zero-range line shapes in the J /ip tt^tt~ and J /ip 7r"'"7r~7r° channels 
had been considered previously by Voloshin [24], but his results were incorrect because of 
conceptual errors related to the treatment of isospin symmetry [10]. 

In order to discriminate between the two binding mechanisms identified above, it is 
necessary to use line shapes that allow for either possibility. There are two criteria that we 
use as guiding principles in constructing the line shapes: 

• unitarity. The line shapes should correspond to multichannel scattering amplitudes 
that are exactly unitary for real values of the interaction parameters. Complex defor- 
mations of the parameters can then be used to take into account effects of additional 
channels that are not treated explicitly. 

• renormalizability. The scattering amplitudes should be derivable from a renormaliz- 
able local quantum field theory. This guarantees that sensitivity to physics at much 
higher energies can be absorbed into the interaction parameters. 

In a local quantum field theory, the production of particles by a short-distance process can 



be represented by local operators. An advantage of a renormalizable field theory is that one 
can construct renormalized operators whose matrix elements are insensitive to physics at 
much higher energies. The line shapes from that production process are then determined 
by the interaction parameters and by the short- distance coefficients of those renormalized 
operators. 

In this paper, we present a solution to the coupled-channel problem for pairs of neu- 
tral charm mesons D*^D^ and D^D*^ and pairs of charged charm mesons D*~^D~ and 
£)+£)*- ^iiat scatter through zero-range interactions and also through an isospin-0 resonance. 
The resulting scattering amplitudes satisfy the conditions of unitarity and renormalizability. 
They depend on 4 interaction parameters and reduce to the Flatte scattering amplitudes of 
Ref. [20] and to the zero-range scattering amplitudes of Ref. [10] in the appropriate limits. 
We also solve the renormalization problem for the local operators that represent the produc- 
tion at short distances of pairs of neutral charm mesons, pairs of charged charm mesons, and 
the resonance. The line shapes for X(3872) produced by the decays B ^ K + X depend on 
the interaction parameters and on the short- distance coefficients associated with the B -^ K 
transition. We assume that the resonance is the Xci(2-P) and we include constraints on the 
interaction parameters and on the B ^ K short- distance coefficients from charmonium phe- 
nomenology. We use our line shapes to analyze data on the X(3872) and data from the Belle 
and Babar collaborations on the production of D*^D^ above the threshold up to 4000 MeV. 
We try to determine whether the data can discriminate between the two mechanisms for 
generating the X(3872) that were described above. 

We begin in Section II by establishing our notation. In Section III, we summarize the 
zero-range and Flatte scattering amplitudes that have been used in previous analyses and 
we present more general scattering amplitudes that include both as special cases. In Sec- 
tion IV, we present experimental and phenomenological constraints on the interaction pa- 
rameters in the low-energy scattering amplitudes. In Section V, we express the line shapes 
of the X(3872) from B meson decays in terms of short- distance coefficients associated with 
the B —^ K transition. In Section VI, we present experimental and phenomenological con- 
straints on the short- distance coefficients. In Section VII, we analyze data on X(3872) and 
data on B ^ K + D*^D^ to see whether they can discriminate between the two binding 
mechanisms. We summarize our results in Section VIII. In an Appendix, we present the 
quantum field theory formulation of the problem of charm mesons that scatter through both 
zero-range interactions and coupling to a resonance. We determine the renormalization of 
the interaction parameters in the scattering amplitudes. We also solve the renormalization 
problem for the local operators that create the charm meson pairs and the resonance. 

II. NOTATION 

The standard isospin multiplets for the charm mesons are {—D~^,D^), {D'^,D~), 
(— D*+,Z}*°), and (Z)*°,D*~), where the first and second states are the upper and lower 
components of the multiplet, respectively. The most natural charge-conjugation phase con- 
ventions are CD° = +D^ and CD*^ = —D*^. The D*D channels with charge conjugation 
quantum number C = + are then [25] 

(D*°D°)+ = ^ (D*°5° - Z}°D*°) , (2a) 

{D*+D-)+ = -^ {D*+D- - D+D*-) . (2b) 



We will refer to {D*^D^)^ and {D*~^D~)^ as the neutral and charged channels, respectively. 
We use concise notation for the masses and widths of the charm mesons. We denote the 
masses of D^, D^, D*°, and D*^ by Mq, Mi, M*o, and M^i, respectively. (The numerical 
subscript is the absolute value of the electric charge of the meson.) The reduced mass 
for D*^D^ is /i = 966.5 MeV. The reduced mass for D*^D^ is larger by about 0.3%, but 
we will ignore this difference. We measure energies relative to the D*^D^ threshold at 
M*o + Mo = 3871.9 MeV. The energy splitting between the D^+D" and D*°Z)° thresholds 
is 

z/n = (M,i + Ml) - (M,o + Mo) = 8.1 MeV. (3) 

The corresponding momentum scale is 



Kii = v^2^z/n = 125 MeV. (4) 

The total width of the D*'^ is measured and the total width of the D*^ can be predicted 
from measurements of D* decays and isospin symmetry [9]. We denote these widths by r^,i 
and r*o, respectively. The PDG value for r*i [26] and the predicted value for r*o are 

r,i = 96 ± 22 keV, (5a) 

r,o = 66 ± 15 keV. (5b) 

The effects of decays of the D*^ and D*^ on charm meson scattering can be partially taken 
into account through energy- dependent widths. In Ref. [9], energy- dependent widths T^q{E) 
and T^i{E) that depend on the energy E of the pair of charm mesons in their center-of-mass 
frame were defined by scaling the physical partial widths for the decays D* — )■ D-k. The 
use of these energy-dependent widths above the D*D thresholds was an error. If the energy 
is above the D*^D^ threshold at -E = 0, the appropriate width for the D*^ is the physical 
width r*o- If the energy is above the D*^D^ threshold aX E = uu, the appropriate width for 
the D*^ is the physical width r*i. Below the D*D thresholds, the energy-dependent widths 
of Ref. [9] give thresholds at the correct energies for the 3-body DDix states. However, 
as pointed out in [17], they do not give the correct energy dependence just above these 
thresholds. Moreover, the effects of the energy-dependent widths on the line shape of the 
X(3872) are numerically small. We will therefore ignore any energy dependence of the 
widths in this paper. 



III. LOW-ENERGY D*D SCATTERING 

In this section, we discuss the low-energy scattering of the charm mesons D* and D. 
We first summarize the universal scattering amplitude of Ref. [9] , which takes into account 
the large scattering length in the neutral channel. We then describe the coupled- channel 
scattering amplitudes of Ref. [10], which take into account zero-range scattering in the 
neutral and charged channels, and the Flatte scattering amplitudes introduced in Ref. [20], 
which take into account scattering through a resonance. Finally, we present more general 
scattering amplitudes that allow for scattering both through a resonance and through zero- 
range interactions. 



A. Universal scattering amplitude 

We first consider neutral charm meson pairs with scattering only in the channel (1^*°/)*')+ 
defined in Eq. (2a). The transition amplitude A{E) for the scattering of nonrelativistically 
normalized charm meson pairs can be written in the form 

A{E) = -f{E), (6) 

where f{E) is the conventional nonrelativistic scattering amplitude expressed as a function 
of the total energy E of the charm mesons in the center-of-mass frame. We measure E 
relative to the D*^D^ threshold. The universal scattering amplitude for an S-wave threshold 
resonance is 

where k{E) = {—2jiE — ieY^'^ and 7 is the inverse scattering length. If 7 is a real parameter, 
f{E) satisfies the constraints of unitarity for a single-channel system exactly. This universal 
scattering amplitude can be derived from a renormalizable nonrelativistic quantum field 
theory with a contact interaction in a single scattering channel. 

The imaginary part of f{E) can be interpreted as a spectral function for the resonance. 
The spectral function associated with the scattering amplitude in Eq. (7) is 

lmf{E) = |/(^)P [Im7 - Im/t(E)] . (8) 

If 7 is a real parameter, the spectral function reduces to 

lraf{E) = 9iE)\fiE)\'./2ilE + 9{^)^5iE + ^yi2fx)). (9) 

There is a threshold enhancement associated with production of D*^D^ and D^D*^ with a 
peak at -E = +7^/(2/i). If 7 > 0, there is also a delta function contribution at E = — 7^/(2/x) 
associated with a bound state with binding energy 7^/(2yu). 

Scattering in the (D*°Z)°)+ channel cannot be exactly unitary, because the D*^ has a 
nonzero width and because the charm meson pair has inelastic scattering channels. Following 
Ref. [9], the dominant effects of decays of the D*^ can be taken into account by replacing 
<E) by 



k{E) = ^-2fiE - tfiT^o, (10) 

where F^q is the width of the D*^ given in Eq. (5b). Following Ref. [9], the effects of inelastic 
scattering channels for the charm-meson pair other than D^D^n^ and D^D^'j can be taken 
into account by making 7 a complex parameter with a positive imaginary part. In the 
expression for the imaginary part of f{E) in Eq. (8), the term proportional to lm.K,{E) is 
the contribution from channels whose ultimate final states are {D^D^7i^,D^D^'y), including 
D*^D^ and D^D*^. The term proportional to Im7 is the contribution from all other channels, 
including J/ipir^TT" . The lm.K{E) term has a threshold enhancement just above the D*^D^ 
threshold. If Re(7) > 0, both terms have a resonant peak just below the threshold that can 
be identified with the X(3872). Its position and width are determined by the pole of the 
scattering amplitude f{E) in Eq. (7). The complex energy of the pole can be expressed as 

^poie = -7V(2/w) - ^r,o/2. (11) 

Quantitative constraints on the real and imaginary parts of 7 are presented in Section IV A. 



B. Zero-Range scattering amplitudes 

In Ref. [10], the universal scattering amplitude of Section III A was generalized to the 
case of coupled channels (D*°Z)°)_|_ and {D*~^D~)^ defined by Eqs. (2) that scatter through 
zero-range interactions. 

1. General case 



We label the two channels {D*'^D^)^ and {D*^D )_|_ by the indices and 1, respectively. 
The transition amplitudes Aij{E) among these two channels define scattering amplitudes 



fij{E) that depend on the total energy E relative to the D*^D^ threshold 



27r 
A,{E) = -UE). (12) 

A'' 

If the charm mesons scatter through a zero-range interaction, the inverse of the 2x2 matrix 
of scattering amplitudes has the form 

i\E)-' = -K-' + K{E), (13) 

where A is 2 x 2 symmetric matrix. The dependence on the energy E is in the diagonal 
matrix 

whose diagonal entries are n{E) = {—2fiE — ieY^'^ and ki{E) = {—2ji{E — uu) — ieY^"^. 
If the three interaction parameters Aqo, Aqi, and An are all real valued, the scattering 
amplitudes fij{E) satisfy the constraints of unitarity for this two-channel system exactly. 
They can be derived from a renormalizable nonrelativistic quantum field theory with zero- 
range interactions. 



2. Isospin symmetry 

The approximate isospin symmetry of QCD reduces the three interaction parameters Aij 
to two independent parameters 70 and 71: 

A = (I/70) Vo + (1/71) Vu (15) 

where the matrices Vq and Vi are projectors onto the isospin and 1 channels: 

n = i ( _; -; ) . (16a) 

p. = i(;;). (16b) 



The scattering amplitudes obtained by inverting Eq. (13) reduce to 

r .p^ -(70 + 7l)+2/ti(E) 

UE) = ^-^^ , (17a) 

UE) ^ ^, (17b) 

_ -(7o + 7i) + 2/s:(i^) 

^''^^^ ME) ' ^^^'^ 

where the denominator is 

DoiE) = 27i7o - (71 + 7o)K(^) + <E)] + 2k,{E)k{E). (18) 

At energies E far from the charm meson thresholds at and uu, the difference between 
K,i{E) and k{E) can be neglected and the isospin symmetry becomes exact. By considering 
this limit, we can identify 70 and 71 as the isoscalar and isovector inverse scattering lengths, 
respectively. We will refer to the model defined by the scattering amplitudes in Eqs. (17) as 
the Zero-Range model. 

3. Optical theorem 

Scattering in the {D*^D^)j^ and {D*^D^)^ channels cannot be exactly unitary, because 
the D*^ and D*^ have nonzero widths and because the charm meson pairs have inelastic 
scattering channels. Following Ref. [10], the dominant effects of decays of D*^ and D*^ can 
be taken into account by replacing k,{E) by the expression in Eq. (10) and by replacing 

K^{E) by 

/ti(^) = V-2/i(E - z/n) - 2/ir,i, (19) 

where r^,i is the width of the D*^ given in Eq. (5a). The contribution of r^,i to the imaginary 
part of Ki{E) is only important near the D*^D~ threshold. Near the D*^D^ threshold, k,i{E) 
is well approximated by the real quantity ku = 125 MeV. Following Ref. [10], the effects of 
inelastic scattering channels other than DDn and DD'-f can be taken into account by taking 
7o and 71 to be complex parameters with positive imaginary parts. 

The imaginary parts of the scattering amplitudes in Eq. (17) can be expressed in forms 
that are consistent with the Cutkosky cutting rules: 

Im^(E) = Im7o Y.ME) Vo,u /r,(i?) + Im7i ^^/.^(i?) Pi,m /;,(i?) 

k,l k,l 

-ImKiE) UE) n^{E) - \mK,{E) fa{E) fl^{E). (20) 

The terms in Eq. (20) proportional to \v[vk{E) and \m.Ki{E) are the contributions from 
channels whose ultimate final states are {DDn, DD'-f), including D*D and DD* . The terms 
proportional to Im7o and Im7i correspond to other inelastic scattering channels with isospin 
and 1, respectively. 

4-. D*^D^ threshold region 

The interaction parameters 70 and 71 can be tuned so that there is a bound state just be- 
low the D*^D^ threshold that can be identified with the X(3872). The scattering amplitudes 

9 



fij{E) have a pole at a complex energy -Epoic that satisfies Eq. (11) with 7 = K^Epoic)- If 
the tiny difference between Ki{Epo\e) and ku is neglected, the vanishing of the denominator 
Do{E) in Eq. (18) reduces to a linear equation for 7, whose solution is 

^ 27170- (71 + 70) /til .21) 

7i + 7o - 2/tii 

An inverse scattering length 7 that is small compared to k,u requires a fine-tuning of 70 and 71 
so that 27170 ~ (71 +7o)/tii. If I71I ^ Kii, 7o must be fine-tuned to near kii/2 ^ +63 MeV. 
Since line shapes near the D*^D^ threshold are extremely sensitive to 7, it is advantageous 
to take 7 to be one of the independent interaction parameters. This can be accomplished 
by eliminating 70 in favor of 7 using 

71^11 + (71 - 2/tii)7 
7o = — 77. ^ , (22) 

(271 - Kii) - 7 

which follows from Eq. (21). 

The scattering amplitudes fij{E) in Eqs. (17) all have poles in the energy variable k{E) 
at k{E) = 7. The residue of the pole of fij{E) has the form Z^ Z- . The ratio of the 
residue factors is 

,71/2 

% = -^^^. (23) 

If the tiny difference between ni (-Epoie) and ku is neglected, the residue Zq is 



1+ '"''"i^' )". (24) 

(71 - f^u) f^iij 

The behavior of the elastic scattering amplitude foo{E) in the entire D*^D^ threshold region 
defined by \E\ <^ uu is dominated by the pole at n{E) = 7. It reduces in this region to 
ZoflE), where f{E) is the universal elastic scattering amplitude in Eq. (7). As 7 — > 0, Zq 
approaches 1. 

C. Flatte scattering amplitudes 

In Ref. [20], Hanhart, Nefediev, and Kalashnikova proposed Flatte line shapes for the 
X(3872) resonance. The Flatte scattering amplitudes can be derived by assuming that 
the coupled channels (D*°Z)'^)+ and {D*~^D~)-^- scatter only through their couplings to a 
resonance. 



1. Isospin symmetry 

If two channels scatter only through their couplings to an isospin-0 resonance, the entries 
of the 2x2 matrix of scattering amplitude defined by Eq. (12) are 

ME) = /11(E) = -foi{E) = /Flattc(i5), (25) 

10 



where fFiatte{E) is the Flatte scattering amphtude: 

^^^^"•^^^^ ^ E - u - {gy2)[K,iE) + KiE)]- ^^^^ 

The threshold functions k,{E) and k,i{E) are given in Eqs. (10) and (19). If we take g to be 
real, the spectral function associated with the Flatte scattering amplitude in Eq. (26) is 



Im/Flatte(^) = |/Flatte(^)P 



2 
-— Imz/ — Im.Ki{E) — lm.K{E) 



(27) 



The notation of Ref. [20] can be obtained by the substitutions 



g^ -^ g, (28a) 

z/ — > Ef-iT{E)/2. (28b) 

In Ref. [20], the effects of the widths of the D*^ and D*'^ were not taken into account. Thus 
the expressions for k,{E) and k,i{E) were Eq. (10) with F^q = and Eq. (19) with F^,i = 0. 
The authors of Ref. [20] did however allow for energy dependence in Imz/, as indicated by 
Eq. (28b). 

2. D*^D^ threshold region 

The interaction parameters z/ and g can be tuned so that there is a bound state just below 
the D*^D^ threshold that can be identified with the X(3872). The scattering amplitude will 
have a pole at a complex energy -Epoie that satisfies Eq. (11) with 7 = /t(_Epoic). If the 
tiny difference between Ki{Epo\e) and ku is neglected, the vanishing of the denominator in 
Eq. (26) reduces to a quadratic equation for 7, one of whose solutions is 



7 = -\g^^i + sJ{\9^^iY - /i(2z/ + g^K^^ + 2F,o). (29) 

The residues of the poles of the scattering amplitudes fij{E) at k{E) = 7 are Z^ Z, , 
where 

zl/' = -Zl/'=(l + ^ + ^Y'\ (30) 

An inverse scattering length 7 that is small compared to ku can be obtained by a fine-tuning 
of z/ such that |z/ + g^K,ii/2\ <C g'^fj,/8,g^Kii/2. The solution for 7 in Eq. (29) then reduces 
to 7 ~ —{2h' + g^Kii +iT^,o)/g'^. In this case, the elastic scattering amplitude foo{E) reduces 
in the entire D*^D^ threshold region \E\ ^ z/n to Zof[E), where f{E) is the universal 
scattering amplitude f{E) in Eq. (7). 

An inverse scattering length 7 that is small compared to ku can also be obtained by 
a double fine-tuning of g and z/ so that g^fi <^ ku and /i|z/| ^ /t^^^. In this case, the 
denominator in Eq. (26) has two zeroes near the D*^D^ threshold. The energy dependence 
of foo{E) in the D*^D^ threshold region is therefore more complicated than the universal 
scattering amplitude. 
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3. Resonance far above D*D threshold 

Alternatively, the parameters v and g can be chosen so that, instead of having a pole 
just below the D*^D^ threshold, the Flatte scattering amplitude in Eq. (26) has a pole 
at an energy well above the D*^D^ threshold. Such a pole could be associated with the 
charmonium state Xci(2-P). If the difference between Ki{E) and k,{E) is neglected, we get a 
quadratic equation for k. The solution for the complex pole is 



E^ -tTj2^u- g^fi - ig^ ^2^// - g^fi^. (31) 

If u has a small imaginary part, the energy and width of the resonance are approximately 

E^ ^ Rez/ — (7^/i, (32a) 

r^ ^ -2Imz/ + 2/v^2^?I(R^^7^^^VyT^V- (32b) 

In the expression for the width, the term with the square root is the partial width for decays 
into D*D. The term proportional to Imz/ corresponds to all other decay channels. As the 

1/2 

energy E^ is increased by adjusting z/, the width T^ increases as E^/ . This behavior is 
characteristic of an ordinary resonance. Near the resonance, the Flatte scattering amplitude 
in Eq. (26) can be approximated by 

/Flatt.(i^)^(l-^ ,^ f^ . ^ ] fBw{E), (33) 

where /bw(-S) is the Breit-Wigner scattering amplitude: 

D. Zero-Range+Resonance scattering amplitudes 

We now generalize the scattering amplitudes of Sections III B and III C to the case of 
coupled channels {D*'^D^)^ and {D*~^D~)^ that scatter through a resonance as well as 
through zero-range interactions. We will refer to the resonance channel as x- 

1. General case 

The transition amplitudes Aij{E) in Eq. (12) define a 2x2 matrix of scattering amplitudes 
fij{E). The expression in Eq. (13) for the inverse of that 2x2 matrix in the case of zero- 
range scattering can be generalized to one that also takes into account the coupling to a 
resonance: 

f{E)-' = -(^A + G^^G^^ +K{E), (35) 

where A is a 2 x 2 symmetric matrix, G is a 2-component column vector, and K{E) is the 
diagonal matrix in Eq. (14). The propagator for this resonance is 

P{E)=[E-iy + G^{A~K-Y^Gy\ (36) 
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If the six interaction parameters Aqo, Aqi, Ah, Gq, Gi, and u are all real valued and if 
r=^o and r*i are set to in the expressions for Eqs. (10) and (19), the amplitudes fij{E) 
satisfy the constraints of unitarity exactly. The scattering amplitudes and the resonance 
propagator can be derived from a renormalizable nonrelativistic quantum field theory with 
two scattering channels and a resonance that interact only through contact interactions. 
The renormalization of this quantum field theory is described in the Appendix. 



2. Isospin symmetry 

The approximate isospin symmetry of QCD reduces the three interaction parameters in 
the matrix A to two independent parameters 70 and 71 defined by Eq. (15). The assumption 
that the resonance has isospin reduces the two interaction parameters in the column vector 
G to a single parameter g: 

Thus the inverse matrix in Eq. (35) reduces to 

1 \ -1 /I 2 \ -1 



A + G- CM = - + -i^ Po + 7i^i- (38) 

E-u J V7o E-uJ 

The scattering amplitudes obtained by inverting the matrix in Eq. (35) can be obtained 
from those in Eq. (17) by replacing 70 by [I/70 + g'^/{E — i^)]^^'- 

ME) ^ |-7.-ro + 2..(E)l(^-. + ,^„)+,^g ^ j3^^j 

/„(B) = i-n -ME -J^+^sh.)+,'-,S ^ ,3g„ 

, ,p> [-71 - 10 + 2k(E)]{E -v + c/jo) + g^g 

ME) = -^^ , (39c) 

where the denominator is 

D{E) = [2-f,-fo-{li + lo)ME) + K{E)]+2K^{E)K{E)]{E-u + g^-fo) 

+g^j^^[-2j, + n,{E) + K{E)]. (40) 

The resonance propagator in Eq. (36) is 

o. r.^ _ 27170 - (71 + lo)[^i{E) + n{E)] + 2k,{E)k{E) 

^^^> - D{E) • ^^^' 

At energies E far from the charm meson thresholds at and uu, the difference between k,i{E) 
and K,{E) can be neglected and the isospin symmetry becomes exact. The scattering am- 
plitudes in Eqs. (39) and the resonance propagator in Eq. (41) depend on four independent 
interaction parameters: 70, 71, g, and z/. We will refer to the model with the scattering am- 
plitudes in Eqs. (39) and the resonance propagator in Eq. (41) as the Zero-Range+Resonance 
model. 
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3. Optical theorem 

The system consisting of the three channels {D*^D^)^^ {D*^D^)^, and x cannot be 
exactly unitary, because the D*^ and D*'^ have nonzero widths, the charm meson pairs have 
inelastic scattering channels, and x iiiay have decay channels other than D*D and DD*. 
We can take into account the dominant effects of decays of D*^ and D*^ by replacing k,{E) 
and K,i{E) by the expressions in Eqs. (10) and (19), respectively. We can take into account 
the effects of inelastic scattering channels other than DDn and DD'-f by taking 70 and 71 to 
be complex parameters with positive imaginary parts. We can take into account the effects 
of decay channels for x other than D*D and DD* by taking z/ to be a complex parameter 
with a negative imaginary part. We choose g"^ to be a real parameter. The imaginary parts 
of the scattering amplitudes in Eq. (39) can be expressed in forms that are consistent with 
the Cutkosky cutting rules: 

Im/,,(E) = f '^"'^l^' 12 1^70- ip. ^''1°'', lmu)Y,UE)Vo,MmE) 
\\E - u + g^-fo\^ \E - u + g^-fol^ /^ -^ 





\E-iy\^ 
- z/ + ^2^0 ^ 


Im7o 




9'\l.? 


E 


E- 


■Z/ + ^2^oP 


Im' 


yi)_^ UE) 


Vi,M 


mE) 





k,l 

-lmK{E) UE) f*^{E) - lmK,{E) U{E) f*^{E). (42) 

The terms proportional to ImnlE) and IinK,i{E) are the contributions from channels whose 
ultimate final states are [DDtt, DD'-f), including D*D and DD*. The terms proportional to 
Im7o, Im7i, and Imz/ correspond to other inelastic D*D scattering channels with isospin 0, 
other inelastic D*D scattering channels with isospin 1, and decay channels of x with isospin 
0, respectively. 



4- Zero-range limit 

In the limits (7 — )■ or z/ — )■ 00, the resonance decouples and the scattering amplitudes for 
the Zero-Range+Resonance model in Eqs. (39) reduce to those for the Zero-Range model 
in Eqs. (17). The Zero- Range model is a good approximation if \E — v\ ^ fi'^|7o|- This 
energy region includes the D*^D^ threshold if \v\ ^ 5*^170 1, in which case the condition on 
the energy reduces to \E\ <^ \u\. 

The Zero-Range model actually has a larger domain of validity if one allows for renor- 
malization of the parameter 70. In the region \E\ <^ |z/|, |z^ — (7^70!, the scattering amplitudes 
in Eqs. (39) reduce to those in Eqs. (17) with the substitution 70 — > ^()V jiy — g'^lo)- 

5. Flatte limit 

In the limits 70, 7i — )■ c>o, scattering proceeds only through the resonance and the scat- 
tering amplitudes for the Zero-Range-|-Resonance model in Eqs. (39) reduce to the Flatte 
scattering amplitudes in Eqs. (25) and (26). The resonance propagator in Eq. (41) reduces 
in this limit to P{E) ^ {—2/g'^)fY\3.ttc{E). The Flatte model is a good approximation in the 
region |7i| ^ ku if the energy satisfies \E\ <^ |7ip//i and \E — u\ <^ g'^\lo\,g'^\li\- This 
energy region includes the D*^D^ threshold if |z/| ^ S'^ItoIjS'^ItiI- 
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The Flatte model actually has a larger domain of validity if one allows for renormalization 
of the couphng constant g. If |7i| ^ Hu and if the energy satisfies \E\ <^ |7ip///, |i/ — g^'^o\ 
and \E — u\ <C K^?^ — z^/7o)7i|, the scattering amplitudes in Eqs. (39) reduce to those in 
Eqs. (25) with the substitution g^ -^ g^ — v/'^o- This energy region includes the D*^D^ 
threshold if |7i| ^ \'~^qu/{u — 5'^7o)|- 

6. D*^D^ threshold region 

The interaction parameters 70, 71, (7, and v can be tuned so that there is a bound state 
just below the D*^D^ threshold that can be identified with the X(3872). The scattering 
amplitudes will have poles at a complex energy -Epoic that satisfies Eq. (11) with 7 = K,{Epoie)- 
The denominator defined in Eq. (40) must vanish at -Epoic- If the tiny difference between 
K,i{Epoie) and /tn is neglected, the equation D{Epoic) = reduces to a cubic polynomial 
equation for 7, which can be written 

[7i/«ii + (71 - 2/€ii)7] [2/i(z/ - 5(^70) + 7^ + ifiT^o] 

= 7o (271 - Kn - 7) (2/i// + 7^ + il^^*o) ■ (43) 

The general behavior of the three poles has been analyzed in Ref. [27] for the case r*o = 0. 
The approximate fine tuning required to obtain a small inverse scattering length can be 
obtained by setting 7 = and F^o = 0: 

7i/tii(z/ - /70) ^ 7o(27i - f^ii)'^- (44) 

If I71I ^ Kii, this fine-tuning condition reduces to Xj^Q—g^jv ^ 2/ ku. An inverse scattering 
length 7 that is small compared to kh can be obtained by fine-tuning 70 to near /€ii/2 ^ 
-|-63 MeV with |z^|/5'^ S> Kh or by fine-tuning u to near —g'^Kii/2 with |7o| ^ Kh. 

Since the line shapes near the D*^D^ threshold are extremely sensitive to 7, it is ad- 
vantageous to take 7 to be one of the independent interaction parameters. The equation 
D{Epoic) = can be solved for 70 as a function of 71, g, u, and 7. If the tiny difference 
between K,i{Epoic) and Ku is neglected, the solution is 

70 = ( (^^^-^^1)-^ + 9' \ ' _ .45. 

V7i'«ii + (71 - 2«;ii)7 iy-Epo\cJ 

If the equation for 7 in Eq. (43) is expanded to first order in the imaginary parts of all 
the variables, one can solve for the imaginary part of 7: 

Im7 ^ - ( ". I Im7o + - I — ) Im7i - - ( ^\ ) (imz/ + |r*o) • 

(46) 
We have simplified the coefficients by setting 7 = 0. We have also used the fine-tuning 
condition in Eq. (44) to express the coefficients in a manifestly positive form. 

The scattering amplitudes fij{E) in Eqs. (39) have poles at n{E) = 7 with residues 
Z^' Zj' . If the tiny difference between fi;i(ii^poie) and Kn is neglected, the ratio of the 
residue factors is 



ryl/2 

Zi 71-7 



Zq^^ 7i - i^u ' 
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(47) 



This ratio does not depend on the resonance parameters v and g. The residue Zq is 

(71 - 7)^7 , 9\l\^^\\ + (71 - 2/tii)7]27" 



1 + 



(71 - /«ll)^/«ll 2/i(7i - /€ll)2(z/ - Epolc^ 



2 



(48) 



As 7 — )■ 0, this residue approaches 1. If the small value of 7 arises from either the fine 
tuning 7o ~ ^^\\|'^ or the fine tuning v ^ —g'^Kii/2, the scattering amplitudes fij{E) in 
the entire D*^D^ threshold region defined by \E\ <^ vn = 8.1 MeV are dominated by the 
pole at k(-E) = 7. They reduce to the universal scattering amplitude f{E) given in Eq. (7) 
multiplied by residue factors: 

f,,{E)^Z]'^ f{E)Z]'\ (49) 

An inverse scattering length 7 that is small compared to kh can also be obtained by a 
double fine-tuning of 70 and z/ so that they satisfy \u\ <^ g^^u, |7o| <^ nn, and \i^ — g'^'Jol <^ 
|z/|. In this case, the residue Zq in Eq. (48) can be significantly smaller than 1. The energy 
dependence of fij{E) can also differ significantly from that of the universal amplitude in 
Eq. (7). 



7. Resonance far above D*D threshold 

One region of parameter space in which the amplitudes simplify is when the resonance 
parameter u is much larger than the energy scale un = 8.1 MeV associated with isospin 
splitting. In this case, the existence of the X(3872) requires the fine-tuning 70 ~ /tii/2. 
In addition to the X(3872) resonance just below the D*^D^ threshold, there is a second 
resonance x '^gH above the D*^D~ threshold. This second resonance could be identified 
with the P-wave charmonium state Xci- 

We first consider this system at energies in the D*D threshold region \E\ < vu. In this 
region, the scattering amplitudes in Eqs. (39) reduce to those in Eqs. (17) for the Zero-Range 
model with the substitution 70 — )■ 7oz//(z/ — g^'jo). 

We next consider this system at energies in the x resonance region. In this region, 
the difference between k{E) and Ki{E) can be neglected and the resonance propagator in 
Eq. (41) reduces to 

""^^^ " E-u+gyiiho-i/mr ('°) 

It has a pole at a complex energy E^ — iT^/2 near u — g^'jo- The position of the pole can 
be calculated by iterating around this approximate solution. The solution to first order in 

7o/«(^x) is 

2 2 

E^ - tTj2 ^ V - (?^o - i^=P^==. (51) 

If 7o and v have small imaginary parts, the energy and width of the resonance are approxi- 
mately 

E^ ^ Re(zy-/7o), (52a) 

r^ ^ -2Imz/ + 2(?2lm7o + 2^'(Re7o)' ^ (52b) 

v/2/xRe(z/ - 5f2^o) 

16 



The last term in Eq. (52b) is the partial width for decays into D*D and DD*. The terms 
proportional to Imz/ and Im7o correspond to other decay channels of x'ci ^^^ ^o inelastic 
isospin-0 scattering channels for charm meson pairs, respectively. As the energy E^ is 

— — —1/2 

increased by adjusting z/, the partial width into D*D and DD* decreases to as E^^ ' . 
This behavior is characteristic of a Feshbach resonance [28]. 

The scattering amplitudes also simplify near the resonance. Using the approximation 

|7o| ^ \2fiE^\^^'^, they reduce to 

foo{E) ^ ME) ^ -U{E) ^ -^ h^{E), (53) 

where /bw(-^) is the Breit-Wigner amplitude in Eq. (34). 



IV. ESTIMATES OF THE INTERACTION PARAMETERS 

The scattering amplitudes for the Zero-Range+Resonance model in Section III D depend 
on the interaction parameters 7, 71, g, and u. In this section, we analyze the constraints on 
these parameters. 

A. Inverse scattering length 

The complex inverse scattering length 7 can be determined from measurements of the 
position and width of the X(3872) resonance in the J /ipn^Ti' decay channel. The position 
of the X(3872) resonance relative to the D*^D^ threshold is given in Eq. (1). There are only 
upper limits on the width of the X(3872) resonance in the J/ipTi^Ti^ decay channel. The 
upper limit from combining the results of the Belle and Babar Collaborations [1, 12] is 

Vx < 2.2 MeV (90% C.L.). (54) 

A lower bound on the width Tx is the width r^,o of the constituent D*°, which is given in 
Eq. (5b). This contribution to the width of the X(3872) can be identified with the decay 
modes {D^D'^n°,D'^D°-f). 

We denote the real and imaginary parts of the complex inverse scattering length by 7re 
and 7im: 

7 = 7re + ilim- (55) 

An alternative pair of variables that can in principle be measured directly are the peak 
position -Emax of the resonance and its full width at half-maximum Ff^hm- The variables 
-Einax and Ff„hm are functions of 71-0, 7im5 and the D*^ width F^,o- They can be expanded in 
powers of F*o [22] : 



i^max = -^-^F,o + ..., (56a) 

^7rc7im . -p, . /t-rv,\ 

+ F*o + (56b) 



fwhm 



/i 



In the expression for Ffwhm in Eq. (56b), the second term F*o can be identified with the 
partial width for decay into (D°l)°7r°, D^D^'-j) while the first term can be identified with 
the partial width into other decay modes. 
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Reference 


Mass (MeV) 


Partial width 


(MeV) 


9 


BG [32] 


3953 


118 




0.37 


ELQ [33] 


3968 


150 




0.40 


BGS [34] 


3925 


165 




0.46 


ELQ [35] 


3920 


81 




0.35 



TABLE I: Results from coupled-channel potential models for the mass of the Xcii^P), its partial 
width into D*D and DD*, and the coupling constant g inferred from Eqs. (32). 

The result in Eq. (1) can be interpreted as a measurement of -Emax- Keeping only the 
leading term in the expansion for -Emax in Eq- (56a), we obtain a determination of 7re: 

7re = 28+^^ MeV. (57) 

The result in Eq. (54) can be interpreted as an upper limit on Ffwhrn- Keeping only the first 
two terms in the expansion for Ff^hm in Eq. (56b), we obtain an upper limit on the product 

of 7i.c and 7im: 

< 7re7im < (32 MeV)l (58) 

B. Resonance parameters 

The resonance parameters are the energy variable u and the coupling constant g. If 
u is much larger than the energy scale Un = 8.1 MeV of isospin splitting, the Zero- 
Range+Resonance model predicts a narrow resonance x whose energy is well above the 
D*+D- threshold. This resonance could be identified with the P-wave charmonium state 
x'ci = Xci(2-P) or a 1 "'""'" cc tetraquark meson or some other 1 "'""'" meson. If x is identified with 
the x'ci! W6 can take advantage of the well-developed charmonium phenomenology based on 
quark potential models to constrain the resonance parameters. 

Predictions from potential models for the mass of x'ci range from about 3920 MeV to 
about 4010 MeV [29-36]. These predictions are all 50 MeV or more higher than the mass of 
the X(3872). Since the effects of couplings of charmonium states to pairs of charm mesons 
are not well understood, we cannot exclude the possibility that they shift the mass of the 
x'ci down into the D*D threshold region. We will take the real part of z/ to be an adjustable 
parameter. 

We can use potential models to estimate the coupling constant g. Conventional potential 
models contain no information about charm mesons. Coupled-channel potential models 
include additional interactions that couple a charm-quark and antiquark to pairs of charm 
mesons. These models can be used to calculate the partial widths for decays of charmonium 
states into pairs of charm mesons. The partial width of x'ci into D*D and DD* has been 
calculated using the ^Pq model [32, 34] and the CCC model [33, 35]. The results for the mass 
of the x'ci ^nd its partial width into D*D and DD* are given in Table I. Although these 
coupled-channel potential models allow scattering of charm mesons, there is no reason to 
expect the scattering lengths to be much larger than the range of charm meson interactions 
in the absence of the fine-tuning of a charmonium state to the D*^D^ threshold. Thus 
the appropriate limit of the Zero-Range+Resonance model in Sec. Ill D is the Flatte limit 
70, 7i ~^ oo, in which the charm mesons scatter only through their coupling to the resonance. 
The energy E^ and the width F^ of the resonance in this limit are given in Eqs. (32). By 
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fitting these expressions for tlie energy and tlie widtli, we obtain tlie coupling constants 
listed in Table I. The average value is 

g = 0.40, (59) 

and the variations are less than 15%. We will use the value in Eq. (59) in the numerical 
analysis in Section VII. 

The effects of other decay modes of the x'ci besides D*D and DD* can be taken into 
account through the imaginary part of the parameter u. The next most important decay 
modes are expected to be the radiative transitions x'ci ~^ J I ''Pi ^"^^ x'ci ~^ V^(2>S')7. The 
partial widths scale like the cube of the photon energy, so they depend on the mass of the 
x'ci- They have been calculated using quark potential models [32-36]. If the mass of the x'ci 
is close to 3872 MeV, the partial widths are roughly 10 keV for J/ip 7 and roughly 60 keV 
for ilj{2S)'~f. These are small enough that the imaginary part of z/ will not have a dramatic 
effect on the line shapes near the D*^D^ threshold. We will therefore set Imz/ = 0. 



C. Charm meson scattering parameters 

The charm meson scattering parameters are 70 and 71. We first consider the natural 
scale for these parameters. Since the low-energy scattering of charm mesons is dominated 
by pion exchange, the obvious estimate for the range of the interaction is l/m^. The 
corresponding estimate of the natural scale for 70 and 71 is m^ ~ 140 MeV. Suzuki has 
suggested that the range of the interaction could be much larger than l/m^ [37]. The 
denominator of the propagator for the exchange of a tt^ of momentum q between D*^ and 
D'^ is g^ + m^o — (M*o — Mq)'^. There is a near cancellation between the terms m'^o and 
(M*o — MqY, because M^,o — Mq is larger than m^o by only 7.1 MeV. If we define an effective 
mass m by expressing the denominator as g^ + m'^, this effective mass is pure imaginary: 
m'^ = —(44 MeV)^. One might be tempted to take |m| as the natural scale for 70 and 71, but 
this would be incorrect. This can be seen by considering the limit m^o — > M^,o — Mq, which 
implies ?72 — )■ 0. In this limit, the potential from the exchange of vr*^ reduces to a 1/r^ potential 
that couples the S-wave and D-wave components of the wavefunction. Although this is a 
long-range potential, it is not a scale-invariant potential. The coefficient C3 in the potential 
C^/r^ provides a scale. In Ref. [38], this coefficient was expressed as C3 = g"^ / {An f^) , where 
g/f-jT is the coupling constant for the D* — Dtt interaction. The corresponding momentum 
scale is (2yuC3)~^ = 360 MeV. This is the natural scale for 70 and 71 when 171 = 0. 

Phenomenological estimates of 70 and 71 can be obtained from any model that can be 
used to calculate the scattering amplitudes for charm mesons. One such class of models is 
meson potential models, in which the degrees of freedom are mesons and their interactions 
are defined by potentials. The simplest such model is one in which the charm mesons interact 
only through the exchange of pions [38-40] . More elaborate meson potential models include 
the effects of the exchange of other mesons [25, 41] or the exchange of quarks [42]. These 
models typically require an ultraviolet cutoff to regularize singularities in the potentials 
at short distances. The scattering amplitudes can be calculated by solving the Schrodinger 
equation. Another class of models that can be used to calculate the scattering amplitudes are 
meson scattering models, which are defined by scattering parameters. In the simplest such 
models, such as the universal theory described in Section III A and the Zero- Range model 
described in Section IIIB, the only degrees of freedom are charm mesons. Other hadrons 
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can also be included as degrees of freedom in meson scattering models. The Flatte model 
in Section III C and the Zero-Range+Resonance model in Section III D include a resonance 
that can be identified with the charmonium state x'ci- More elaborate meson scattering 
models that include many other mesons as degrees of freedom have also been considered 
[43] . In meson scattering models, the scattering amplitudes are calculated by solving integral 
equations, such as the Lippmann-Schwinger equation. These models typically require an 
ultraviolet cutoff to regularize singularities at large momenta. If the model is renormalizable, 
like those described in Sections III A, IIIB, IIIC, and HID, all dependence on the cutoff can 
be absorbed into the scattering parameters. Although many meson potential models and 
meson scattering models have been applied to the X(3872), charm meson scattering lengths 
have not been calculated in these models. Thus the published results for these models do 
not provide any direct estimates of the scattering parameters 70 and 71. 

A particularly convenient observable for constraining 70 and 71 is the ratio of the residues 
of the poles in the elastic scattering amplitudes foo{E) for D*^D^ and fu{E) for D*^D~ 
at the X(3872) resonance. We denote the residue of the pole in fij{E) at k,{E) = 7 by 

Z^ Zj . The quantity \Zi/Zo\^/'^ can be interpreted as the ratio of the wavefunctions 
at the origin for the {D*+D-)+ and (L'*°D°)+ components of the X(3872) [44]. For the 
Zero-Range+Resonance model, the ratio Z^ /Zq is given in Eq. (47). It is determined 
primarily by 71, because I7I ^ nu- It does not depend on the resonance parameters g and 
u. This ratio has been calculated in a meson scattering model whose degrees of freedom are 
SU{4) multiplets of pseudoscalar and vector mesons [43]. The parameters of the model were 
fine-tuned so that the binding energy of the X(3872) is 0.4 MeV. The absolute value of the 
ratio of the residues was calculated to be \Zii/Zqq\^^'^ = 0.9923. If we insert the expression 
for the ratio of the residues in Eq. (47), we obtain two solutions for 71: +76 MeV and 
— 12, 500 MeV. We expect 71 to be large compared to ku = 125 MeV and positive, because 
the pion-exchange interaction is repulsive in the C = + isospin-1 D*D channel. Neither of 
the two solutions for 71 are consistent with this expectation. It is possible that the model 
of [43] does not describe scattering in the isospin-1 channel with sufficient accuracy. 

V. LINE SHAPES OF X(3872) IN B MESON DECAY 

In this section, we consider the line shapes of the X(3872) resonance in the decay B — )■ 
K + X. We first present a general formulation of the problem of production by a short- 
distance process in terms of an effective field theory that describes the D*D threshold 
region. We summarize the universal line shapes of Ref. [9], which take into account only the 
large scattering length in the neutral channel. We summarize the line shapes of Ref. [10], 
which take into account zero-range scattering in the neutral and charged channels. We then 
describe the Flatte line shapes introduced in Ref. [20], which take into account scattering 
through a resonance. Finally, we present more general line shapes that allow for both 
scattering through a resonance and through zero-range interactions. 
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A. Effective field theory formulation 

Our starting point for the derivation of the hne shapes produced by the decay B^ — )■ 
K'^ + X is the optical theorem for the width of the B^: 

nB^] = ^lmA[B^^B^], (60) 

Mb 

where iA[B^ — )■ B~^] is the one-meson-irreducible forward amphtude for 5+ at leading 
order in the electroweak interactions and to all orders in QCD interactions. The imaginary 
part of this amplitude has a contribution from the intermediate state K^ + X(3872). In 
addition to the X(3872) itself, there are other sets of particles with the quantum numbers 
J^*^ = 1"*"^ that have enhanced production rates near the D*^D^ threshold. We will denote 
these states collectively by the symbol (1^^). We will use the phrase X(3872) resonance to 
refer specifically to the peak in the energy distribution just below the D*^D^ threshold. 

Decays of the 5+ proceed through weak interactions that are mediated by the W boson. 
Because the mass of the W is so much larger than that of the B meson, the decays can 
be described completely within QCD using effective field theory methods. The effects of 
W exchange can be reproduced by an effective weak Hamiltonian. The leading terms that 
contribute to B^ — )> K^ + X are current-current operators: 

"Hweak = Ccc h^il - 75)C C7^(l " 75)s + Cnc h'^i^ " Ts)^ C7^(l " 75)c, (61) 

where Ccc and Cnc are short-distance coefficients. The forward amplitude in Eq. (60) can 
be expressed as the expectation value in the B^ meson of the Fourier transform of a bilocal 
QCD operator: 

A[B+ ^ B+] =1 j d^x (5+|T?/t^^,(x)?/„eak(0)|5+). (62) 

The K^ and X(3872) produced by the decay of B^ have recoil momenta of 1140 MeV, 
which is large compared to the momentum scale associated with the resonance. If the 
X(3872) can be described using an effective field theory (EFT) for charm mesons, the decay 
B^ — )> X+ + X(3872) can also be described completely within that EFT. In the rest frame of 
the X(3872), the B'^ — )■ K^ transition acts like a point source for pairs of charm mesons and 
possibly other degrees of freedom described by the EFT. As far as the resonance is concerned, 
the effects of the B^ — ?> K^ transition can be reproduced by local EFT operators O™" acting 
on the EFT vacuum. The superscript m is a Cartesian vector index and the subscript i 
labels the various operators. The contribution to the forward amplitude in Eq. (62) from 
transitions to K^ that create 1"^"^ states with invariant mass in the D*D threshold region 
can be expressed in terms of expectation values in the EFT vacuum of the Fourier transforms 
of bilocal EFT operators: 



(27r)4 Pl-m\ + ie 



«j 



xi Idt e'^* jdh (0|TC;^(r,t)Cr(O,0)|0). (63) 

The energy E is defined by (Pg — Px)^ = (^*o + Mq + E)'^, where Pb and P^ are the 4- 
momenta of the B~^ and K~^. The integral over the 4-momentum P^ of the K'^ is implicitly 
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restricted to the region near the D*D threshold where E is small. The short-distance 
coefficients (7*'"^ specify the linear combinations of EFT operators ^^ (ji,mQm (^^^\^]^ qj^ 
implied sum on the vector index m) that reproduce the effects of the B^ — )■ K^ transition. 
The short-distance coefficients C**'"^ can in principle be determined by matching matrix 
elements of "Hweak in QCD with matrix elements of ^^ (jt,mQm -^^ ^j^g EFT. If we consider 
a D*^D^ state in the D*D threshold region, the matching condition is 

(ir+Z}*°5°|?/„eak(0)|fi+) = ^C*'™(D*°D°|Cr(O,0)|0). (64) 

i 

Although we cannot calculate the QCD matrix element on the left side, we can use the 
Lorentz invariance of QCD to determine the tensor structure of the coefficients C*'™. We 
choose the charm meson pair D*^D^ to have zero relative momentum. The QCD matrix 
element on the left side of Eq. (64) must be linear in the polarization 4- vector e^ of the D*^ 
and it can depend only on the external 4- momenta Pb-, Pk, and Px = Pb — Pk of the B~^, 
K^, and D*^D^ pair. Since the operator "Hweak is a Lorentz scalar and Px ■ e = 0, the QCD 
matrix element must be Pb ■ e multiplied by a constant. The EFT matrix element on the 
right side of Eq. (64) must be linear in the polarization 3-vector e of the D*^ in the D*^D^ 
center-of-mass frame and there are no 3-momenta that it can depend on. Since the operator 
(9™ is a Cartesian vector with index m, the EFT matrix element on the left side of Eq. (64) 
must be e™" multiplied by a constant. By matching the two sides of Eq. (64), we conclude 
that the short-distance coefficients must have the form 

^i,m ^ _^c^^''pj^L"^ (65) 

where the coefficients C^+ '* are complex numbers and L™ is the boost matrix from a general 

Lorentz frame to the center-of-mass frame of the D*'^D^ pair. The prefactor M^ has been 
inserted for later convenience. The boost matrix satisfies Px^^ = and the tensor identities 

^mn^rnj^n ^ -g^^ + P^^P^J P^, (66a) 

g^'^L'^L^ = -5™". (66b) 

The imaginary part of .4.(1++) can be obtained from Eq. (63) by using cutting rules. The 
relevant cuts run through the K~^ propagator and through the EFT Green's function: 

ImAi++) = ^E^f-^''(C'f:'^r||^27r5(P^-mi)PKL-L:: 






Imi fdte'^' fdh- ^0|TC"^(r,t)Cr(O,0)|0y (67) 

Rotational invariance implies that the EFT matrix element is proportional to (5™". Using 
Eq. (66a), the integrand of the integral over P^ in Eq. (67) becomes Lorentz invariant. It 
can be reduced to an integral over {Pb — PrY or, equivalently, E: 
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In the D*D threshold region, the energy E is neghgible compared to M*o + Mq. The 
amphtude in Eq. (67) therefore reduces to 

where Fij{E) is the EFT matrix element: 

Fij{E) = t fdt e'^* fd^r (0|TO™^(r, t)Cr(0, 0)|0). (70) 

Note that there is an implied sum over the repeated Cartesian vector index m. Inserting 
Eq. (69) into Eq. (60), we get a factorization formula for the inclusive energy distribution 
summed over l"*""*" states: 



^[B^ ^K^ + (1++)] = rf J]Cj:'^(Cj:'^r ImF,,(i?). (71) 

The prefactor has dimensions of energy: 

^K ^ X'/\MB,mK,M,o + Mo) 
^ 967r2M|(M,o + Mo) ' ^ ^ 

If the short-distance coefficients C^+ '* are chosen to be dimensionless, the functions Fij{E) 
have the same dimensions (momentum)"^ as a scattering amplitude. 

B. Universal line shapes 

In Ref. [9], Braaten and Lu derived expressions for the universal line shapes for 1^^ states 
near the D*^D^ threshold produced by the decays B ^ K + (1^^) that take into account 
the large scattering length in the neutral channel {D*^D^)^. The D*^D^ threshold region 
can be described by an effective field theory for the neutral charm mesons D'^, D*^, D^, and 
Z)*° in which the only interaction is zero-range scattering in the {D*^D^)^ channel. The 
scattering amplitude f{E) is given in Eq. (7), with k,{E) given by Eq. (10). It depends on 
a single interaction parameter 7. 

The effects of the B^ — )■ K^ transition can be reproduced in this effective field theory by 
a local operator (9™ with a Cartesian vector index m that creates a pair of charm mesons 
in the (D*°Z)°)+ channel. The Green's function defined by Eq. (70) for this operator is 
proportional to the scattering amplitude f{E). The normalization of the operator Oq can 
be chosen so that Foo{E) is exactly f{E). The resulting factorization formula for the inclusive 
energy distribution of l"*"^ states produced by the transition B^ — )■ K^ is 



— [5+ ^K+ + (1++)] = rf |c5;'Tlm/(E), (73) 



where C^+ ' is a complex short- distance coefficient. The imaginary part of the scattering 
amplitude in Eq. (73) is given by the optical theorem in Eq. (8). The term proportional to 
Im.K{E) is the contribution from channels whose ultimate final states are {D^D^tt^, D^D^'-f). 
The term proportional to Im7 is the contribution from all other channels. 
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The inclusive energy distribution for 1 "'""'" states produced by the transition B^ — )■ K^ is 
given by an expression identical to Eq. (73) except that C^+ ' is replaced by an independent 
short-distance coefficient C^o ' ■ Since the rate in Eq. (73) does not depend on the phase of 

C^+ ' , the short-distance coefficients for the B^ — )> K^ and B^ — )• K^ transitions depend 
on only two independent real constants. 



C. Zero-range line shapes 

In Ref. [10], the universal line shapes in Section VB were generalized to the Zero- Range 
model of Section IIIB. This model can be derived from a renormalizable quantum field the- 
ory for neutral and charged charm mesons with zero-range interactions among the channels 
[D*^D^)^ and {D*^D^)-^.. The effects of the B~^ — )■ K~^ transition can be reproduced by a 
linear combination of two local operators that create pairs of charm mesons in the (D*°Z)°) + 
and (D*+D^)+ channels. The two operators mix under renormalization. For an appropriate 
choice of the renormalized operators Oq and Oi, the Green's functions Fij{E) defined by 
Eq. (70) are simply the scattering amplitudes fij{E). The resulting factorization formula 
for the inclusive energy distribution of l^"*" states produced by the transition i?+ — )■ K~^ is 

^[5+ ^K^ + (1++)] = rf J2CB"'\C^"'n^UE), (74) 

where the sums are over i,j G {0, 1}. The short-distance coefficients C^+ '* can be complex. 

Since the rate in Eq. (74) does not depend on a common phase in C^+ ' and C^+ ' , the short- 
distance coefficients are determined by three independent real parameters. The scattering 
amplitudes fij{E) are given in Eqs. (17) and expressions for their imaginary parts are given 
in Eq. (20). Note that Eq. (20) gives two different decompositions of the imaginary part of 
the function f\o{E) = fm{E) into terms hnear in the imaginary parts of k{E), k,i{E), 70, and 
7i. To obtain an expression for dT/dE that is consistent with the Cutkosky cutting rules, 
these two different expressions must be inserted for Im/oi(-E) and Im/io(-E) in Eq. (74). 

The inclusive energy distribution for 1++ states produced by the transition B^ — ?> K'^ is 
given by an expression identical to Eq. (74) except that the short-distance coefficients C^+ '* 

are replaced by C^o'*- Since the effective weak hamiltonian "Hwcak in Eq. (61) is an isospin 
scalar, the operator ^^ (jt,mQrn j^^^g^ qI^q Iqq q^^i isospin scalar. The approximate isospin 

symmetry of QCD can then be used to relate C^o '* to C^+ '* [10]: 

C^o ' = (~'B+ ' ' (75a) 

Thus the two short-distance coefficients for the B^ — )■ K^ transition are determined by the 
same three independent real constants as the short-distance coefficients for the B~^ — )■ K~^ 
transition. In Ref. [10], the expressions for C^o '* had the opposite signs from those in 
Eqs. (75). This error did not affect any of the physical results in Ref. [10]. 
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D. Flatte line shapes 

In Ref. [20], Hanhart, Nefediev, and Kalashnikova proposed line shapes for the X(3872) 
resonance in the D*^D^ threshold region that correspond to the Flatte parameterization of 
the D*^D^ elastic scattering amplitude in Eq. (26). Their line shapes in the D*^D^ and 
J/ip TT+TT^ decay channels are 

"^^ [B^ ^ K^{D*'D'^,D'D*')] = ^^^^\M^,,,{E)\'y^2ilE 9{E), (76a) 



dE^ /"^ J ^g 



[B+ ^ ir+J/^TT+TT-] = _^., V Flatte(i^)P Tjj^^^{E), (76b) 



where r[i?+] is the total width of the B^ and i3 is a short-distance coefficient associated 
with the B^ — )> K^ transition. The energy dependence of the function Tji^t^^{E) is the 
decay rate of a resonance of energy E into J /iI)'k^ti~ through the decay into J /ip p followed 
by p — )■ 7r+7r~. 

The Flatte line shapes can be derived in the effective field theory framework of Section V A 
by assuming that the short-distance process is the creation of the resonance through an op- 
erator (9™. The Green's function F22{E) is proportional to the propagator of the resonance. 
The constant factor can be chosen so that F22{E) is equal to the Flatte scattering amplitude 
/Fiatte(-E) in Eq. (26). The inclusive energy distribution in Eq. (71) is then 



^[5+ ^K^ + (1++)] = rf |Cj;'Ylm/Mattc(i?), (77) 



where C^+ ' is the short- distance coefficient for the operator that creates the resonance. 
The imaginary part of the Flatte scattering amplitude is given in Eq. (27). 

The line shapes in Eqs. (76a) and (76b) can be obtained by replacing r^\C^^' p by 
2Br[B~^]/{TTg'^) and by making the substitutions 

-Im/t(E) — > ^2pE 9{E), (78a) 

_Inizy — > Tj/^^^{E)/2. (78b) 

The line shape for D*^D^ in Eq. (76a) vanishes for E < 0. Thus the substitution in Eq. (78a) 
discards all contributions to [D^D^-k^^ D^D^'-j) final states from the X(3872) resonance below 
the D*^D^ threshold. 



E. Zero-range-|-resonance line shapes 

We now generalize the line shapes of Sections V C and V D to the Zero-Range+Resonance 
model of Section HID, in which the coupled channels (D*°Z)°)+ and {D*^D^)^ scatter 
through a resonance as well as through zero-range interactions. 



1. General case 

In the Zero-Range+Resonance model, the scattering amplitudes fij{E) are given by 
Eq. (35) and the resonance propagator P{E) is given in Eq. (36). They depend on the 
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three independent entries of the symmetric matrix A, the two components of the column 
vector G, and u. As described in the Appendix, this model can be derived from a renormal- 
izable field theory for the neutral and charged charm mesons and the resonance x- 

The effects of the B^ — > K^ transition can be reproduced in this model by a linear 
combination of three local operators: two of them create pairs of charm mesons in the 
[D*^D^)^ and {D*^D^)^ channels and the third creates the resonance. The solution to 
the renormalization problem for these operators is presented in the Appendix. Two of the 
renormalized operators Oq and O" can be chosen so that the Green's functions Fij{E) 
defined in Eq. (70) for i,j = 0, 1 are just the scattering amplitudes: 

F,,{E) = f,,{E), ^,JG{0,1}. (79) 

The remaining operator O2 can be chosen so that F22{E) is the resonance propagator P{E) 
given in Eq. (36) multiplied by a constant: 

F22iE) = -{G^G/2)P{E). (80) 

The multiplicative factor of G^G has been inserted only to ensure that F22{E) has the same 
dimensions as a scattering amplitude. The two remaining Green's functions Fq2{E) and 
Fi2{E) are then given by 

F,2{E) = (G^G/2)V2 (^(^) [(^ _ ^)^ ^ ^^T]-i ^^^ ^81) 

The inclusive energy distribution for 1^+ states produced by the transition B^ — )■ K~^ is 
given by the factorization formula in Eq. (71), where the sums are over i,j G {0, 1, 2}. The 
short-distance coefficients C^+ '* can be complex. Since the rate in Eq. (71) is not affected 

by a common phase in C^+ '*, these three coefficients are determined by 5 independent real 
constants. 



2. Isospin symmetry 

Using the approximate isospin symmetry of QCD, the matrix A can be expressed in terms 
of two independent parameters 70 and 71 using Eq. (15) and the column vector G can be 
expressed in terms of a single coupling constant g using Eq. (37). The scattering amplitudes 
fij{E) are given in Eqs. (39). The remaining Green functions Fij{E) in Eqs. (81) and (80) 
reduce to 

F„,(E) = "'^'°'7(+;''^" . (82a) 

F.(E) = -£!2£hll±^!M. (82b) 

where the denominator D{E) is given in Eq. (40) and Dq{E) is given in Eq. (18). 

The inclusive energy distribution for 1"*"^ states produced by the transition B^ — )> K^ is 
given by a factorization formula identical to Eq. (71) except that the short- distance coeffi- 
cients C^+ '* are replaced by G^o '*• Isospin symmetry implies the relations between G^o '* 
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and C^+ '* for i = 0, 1 in Eqs. (75). The assumption that the resonance has isospin imphes 
a similar relation between the short-distance coefficients for i = 2: 



C 



K+,2 
B+ 



c 



K0,2 
BO ■ 



(83) 



Thus the three complex short- distance coefficients for the B^ — )■ K^ transition are deter- 
mined by the same five independent real constants as the short-distance coefficients for the 
B^ — )> K^ transition. 



3. Optical theorem 

The inclusive energy distribution in Eq. (71) depends on the imaginary parts of the 
amplitudes Fij{E). Expressions for the imaginary parts of Fij{E) for i,j G {0, 1} that are 
consistent with the Cutkosky cutting rules are given in Eqs. (42). We also need expressions 
for the imaginary parts of Fq2{E), Fi2{E), and F22{E) that are consistent with the Cutkosky 
cutting rules. The imaginary parts of Fi2{E) = F2i{E) can be expressed as 



lr^Fi2{E) 



-\m[Uo{E) - hi{E)] 



7o 



+ '-[UE)-UE)] 



■"-^•(-' - y U^i^ 



+ 



E-v 



E -u + g'^'jo 

E-u 
E - u + g'^jol^ 

lm[UE)-ME)] 



lmjo + 



7o 



E - iy + g^-fol^ 



Imz/ 



(84a) 



9' 



2 \\E-iy + g^^o\ 



-Im7o + 



7o 



i£_„+,.-,^H'^«<^'-^"''^" 



(84b) 



The imaginary parts of fiQ^E), fn^E), foi{E), and fii{E) can be decomposed into terms 
linear in the imaginary parts of k,{E), k,i{E), 70, 71, and u by using Eqs. (42). Note that 
Eqs. (84a) and (84b) give two different decompositions of the imaginary parts of the function 
Fi2{E) = F2i{E). To obtain an expression for dV/dE that is consistent with the Cutkosky 
cutting rules, these two different expressions must be inserted for Im.Fi2{E) and lTaF2i{E) 
in Eq. (71). The imaginary part of F22{E) can be expressed as 



ImF22(^) = 2\F, 
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-—\m.v 
9 



-7i[/t(E) + Ki{E)] + 2ki{E)k{E) 



Do{E) 



Im7o 



^oHE)-miE)] 



-2 



7o[-7i + /«i(^)] 



Im/t(E) - 2 



DoiE) 
where Dq{E) is the denominator given in Eq. (18). 



70 [-71 + <E)] 



DoiE) 



Do{E) 
lm.K,i[E) 



Im7i 

(85) 
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4-. D*^D^ threshold region 

The interaction parameters 70, 71, ly, and g can be tuned so that there is a bound state 
just below the D*^D^ threshold that can be identified with the X(3872). The amplitudes 
Fij{E) have poles at k,{E) = 7 with residues Z/ Zj . The residue factors Zq and Z/ 
are given by Eqs. (47) and (48). The ratio of Z2 and Zq is 

zl^^ _ g^bitiu + (71 - 2/€ii)7] .ggv 



Zy2 2(7i-«:ii)(z/-i5, 



^0 



pole J 



If the small value of 7 is obtained by a single fine-tuning of 70 or u, the behavior of these 
amplitudes in the entire D*^D^ threshold region \E\ ^ vn ^ 8.1 MeV is dominated by the 
pole at k{E) = 7. They all reduce to the universal scattering amplitude f{E) in Eq. (7) 
multiplied by residue factors: 

F,,{E) ^ Z\l^ fiE) Z]l\ (87) 

The inclusive energy distribution in Eq. (71) simplifies in the D*^D^ threshold region: 

dV 



K 



— [5+^ir+ + (i++)]^r^ 



^K+fi _ 7i nK+,1 _ 9 71^11 r^K+,2 
^' 7i-«:ii ^" 2(71 -/.n)^^" 



2 
ImfiE), 



where lraf{E) is given in Eq. (8). The coefficients Z, of C^+ '"' have been simplified by 
taking the limit 7 — )■ and by neglecting the imaginary parts of the interaction parameters. 
The result in Eq. (88) agrees with the universal factorization formula in Eq. (73) after 
a redefinition of the short- distance coefficient C^+ ' . By using the expression for Im7 in 
Eq. (46), the energy-dependent factor in Eq. (88) can be expressed as a linear combination 
of the imaginary parts of 70, 71, i^, and k,{E). There are additional contributions to the 
coefficients of the imaginary parts of 70, 71, and u that come from cuts through the factors 
of Z^' and Zj' in Eq. (87), but they are suppressed compared to those in Eq. (88) by 
factors that are at least as small as 7//tii or k{E)/ Kh. 



5. Resonance far above D*D threshold 

In the region of parameter space in which u is much larger than the energy scale un = 
8.1 MeV of isospin splitting, there is a second resonance x well above the D*^D~ threshold 
in addition to the X(3872) resonance just below the D*^D^ threshold. In this case, the 
amplitudes Fij[E) can be simplified. 

For energies E in the x resonance region, the line shapes are dominated by a Breit-Wigner 
resonance whose energy and width are given by Eqs. (52). The scattering amplitudes fij{E) 
reduce near the resonance to the expressions in Eq. (53). The remaining amplitudes reduce 
to 

F,,iE) ^ -FuiE) ^ -^P^ fBw{E), (89a) 

F22{E) ^ /bw(^), (89b) 
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where fswiE) is the Breit-Wigner scattering amphtude in Eq. (34). Note that Fij{E) 
differs from /bw(-E) by a factor of —i'-yo/ ^y2fiE^ for every subscript and by a factor of 
+i'yo/ ^,/2jiE^ for every subscript 1. Thus the expression in Eq. (71) for the inclusive energy 
distribution reduces near the resonance to 



dT 

dE 



[5+ ^K+ + (1++)] ^ r 
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yX+,0 



K+,1-. 



,/2ilE, 
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(90) 



Upon integrating over the energy, the partial width from the Breit-Wigner resonance is 

2 



r[5+ -^ K^ 



Xcl\ 






V^i^: 
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-c-,-) + c^,. 



(91) 



6. Zero-range limit 

If the parameters satisfy |z/| ^ S'^ItoI, the Zero-Range+Resonance model reduces to the 
Zero- Range model in the energy region \E\ <^ |z/|. The scattering amplitudes reduce to 
those in Eqs. (17). The amplitudes Fq2{E), F^i^E), and F22{E) reduce to 



Fq2{E) 

Fu{E) 

F22{.E) 



lyDoiE) 

g^loi-li + t^{E)] 



uDo{E) 



2u' 



(92a) 
(92b) 
(92c) 



They are suppressed by a factor of g^'yo/u. Thus the hne shapes reduce to those of the 
Zero- Range model in Eq. (74). 

7. Flatte limit 



If the parameters satisfy |7i| 3> ku and z/ <^ S'^ItoIjA'^ItiI, the Zero-Range+Resonance 
model reduces to the Flatte model in the energy region \E\ <^ llil'^ / 1^, g'^llol, g'^llil- The 
scattering amplitudes fij{E) reduce to those of the Flatte model in Eqs. (25). The remaining 
amplitudes Fo2{E), Fi2{E), and F22{E) reduce to 



Fo2{E) ^ -Fu{E) ^ F22{E) ^ hi^E)- 



(93) 



where fF\a,tte{E) is the Flatte scattering amplitude in Eq. (26). Note that Fij{E) differs from 
/Fiattc(-E') by a factor of —1 for every subscript 1. Thus the expression in Eq. (71) for the 
inclusive energy distribution reduces to 



dr_ 

dE 



[B+ -^ K^ 



-)] 



^ B 



^B+ ~ ^B+ + ^B+ 



Im/Flatte(-E). 



(94) 
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The imaginary part of the Flatte scattering amphtude is given in Eq. (27). The inclusive 
energy distribution in Eq. (94) agrees with that for the Flatte model in Eq. (74) after a 



redefinition of the short-distance coefficient C^+ ' . 



If there is a resonance well above the D*^D^ threshold that can be identified with x'ci^ 
the Flatte scattering amplitude fpiatteiE) in Eq. (26) reduces near the resonance to the 
Breit-Wigner amplitude in Eq. (34). Upon integrating over the energy, the partial width 
from the Breit-Wigner resonance is 



ng 



2 



K+,0 ^K+,1 , ^K^ 



T[B^ ~,K^ + x^J ^ '^V- Cll '^ - Cli '^ + C^; '^ . (95) 



VI. CONSTRAINTS ON THE B ^ K TRANSITION COEFFICIENTS 

In this section, we analyze the constraints on the short- distance coefficients for the B ^ K 
transitions from information about B meson decays. 

A. B meson decays into K -\- X(3872) 

Products of the branching fractions for the decay B -^ K + X(3872) followed by the 
decay of X(3872) into specific final states have been measured by the Belle and Babar 
Collaborations. The most precise measurements are for the final state J/ipir'^TT" [12, 13] 
and for the final states D°D°7r°_ [7] and D°D°7 [15, 16]. We denote the final states D°D°7r° 
and D^D^'-j collectively by [D^D^ti^, D^D^'-y). The measurements for the final state D^D^tx^ 
in Ref. [7] can be interpreted as measurements for {D^D'^TT^,D^D^'~f) multiplied by a 62% 
branching fraction. Since D*'^ ultimately decays into D°7r° or -0^7, the measurements for 
the final states D*^D^ and D^D*^ in Refs. [15, 16] can be interpreted as measurements for 
{D'^D^7r^,D^D^'y). The measurements for 5^ decays are more precise than those for B'^ 
decays. We introduce a simple notation for the product of branching fractions in B~^ decays: 

Br+ [final] = Br[fi+ ^ K+ + X(3872) ^ K+ + final]. (96) 

Another convenient observable is the ratio Ro+ of the products of branching fractions for 
B^ decays and B~^ decays: 

,^ _ Br[go ^K' + X(3872) ^ K^ + final] 
"^'^ [^""^] = Br[i?+^;^+ + X(3872)^ii-^ + final] - ^^^^ 

The measurements of Br+ and -Ro+ are summarized in Table II. To obtain the average 
values of the measurements for J/ipn^Tv^ and {D^D^ti^, D^D^'-y) in Table II, statistical and 
systematic errors were added in quadrature. The ratio of i?o+ for {D^D^ir, D^D^'j) and i?o+ 
for J /ipTi^TT" is 2.1 ± 0.9, which differs from 1 by more than 1 standard deviation. The 
ratio of Br+ for {D^D^tc, D^D^'-f) and Br+ for J/ipTT^ii^ is 11.0 ± 2.5. For reasons that will 
be discussed below, this result does not necessarily imply that the branching fraction for 
X(3872) into {D^D^TT,D^D^'~f) is an order of magnitude larger than for J/iI)ti^ti~ . 

In the Zero-Range+Resonance model, the products of branching fractions can be cal- 
culated in terms of the scattering parameters 7, 71, g, and v and the B -^ K transition 
coefficients C^+ '*, i = 1, 2, 3. The inclusive energy distribution for 1^+ states is given in 
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Reference 


Decay mode 


Br+ X 10^ 


^0+ 


Babar [12] 


J/'ll>'K^TT^ 


8.4 ±1.5 ±0.7 


0.41 ± 0.24 ± 0.05 


Belle [13] 


J/lfj-K^TT^ 


8.10 ± 0.92 ± 0.66 0.82 ± 0.22 ± 0.05 


Average 


J/lfjlT^'K^ 


8.2 ±0.9 


0.63 ±0.17 


Belle [7] 


D^D^'^ 


102 ± Sltll 


1.63 ±1.03 


Babar [15] 


D*^D^,D^D*^ 


167 ± 36 ± 47 


1.33 ±0.69 ±0.43 


Belle [16] 


D*^D°,D^D*° 


77 ± 16 ± 10 


1.26 ±0.65 ±0.06 


Average 


{D^D^7r^,D^D^^) 


90 ±17 


1.35 ±0.46 



TABLE II: Data on products of branching fractions for B — )■ i^ + X(3872) followed by the decay of 
X into specified decay modes. The first column gives the reference for the experimental result by 
the Belle or Babar Collaboration. The second column specifies the decay mode for X(3872). The 
third column gives the product Br-|_ of the branching fractions for B~^ decays. The fourth column 
gives the ratio i?o+ of the products of the branching fractions for B^ decays and B~^ decays. Results 
with both a statistical and a systematic error are experimental results. Results with a single error 
were obtained by combining statistical and systematic errors in quadrature. 



Eq. (71). The imaginary parts of the amplitudes Fij(E) are given in Eqs. (42), (84), and 
(85). The contribution from [D^D^ti^ ,D^D^'-)) is the sum of all the terms with the factor 
\m.K{E). Since the decay of X (3872) into J /ip 7i~^7r~ is dominated by the decay into J/ip and 
a virtual p", this decay mode has isospin 1. The contributions from all isospin-1 final states 
other than DDn and DD'-f is the sum of all terms with the factor Im7i. The contribution 
from J/il)Ti^Ti~ is obtained by replacing Im7i by a term (Im7i)'^/'^'^ ^ . 

We will use several results from Table II to constrain the B ^ K transition coefficients. 
We will use the average values of the ratios -Ro+: 

i?o+[J/^7r+7r"] = 0.63 ±0.17, (98a) 

i?o+[/^°5V,D°D°7] = 1.35 ±0.46. (98b) 

We will also use the product of the branching fractions for B^ to decay into K^ and D^D^n^ 
or D^D^-f. 

Br+[D°i)°7r°, D°i)°7] = (9.0 ± 1.7) x 10"^ (99) 

We will not use the product of the branching fractions for B^ to decay into K~^ and 
J/^7r~^7r~, because it would introduce the additional unknown parameter (Im7i)"'/^'^ ^ . 

We ffist consider the constraints on the parameters from i?o+ [J/'ip Tr"*"?!"] . Since J/ip 7r"'"7r~ 
has isospin 1, the only resonant enhancement comes from the X(3872) resonance just below 
the D*^D^ threshold. The dominant contribution to the product of branching fractions will 
come from the D*^D^ threshold region, where the line shape from the decay of B^ can be 
approximated by the simple expression in Eq. (88). The line shape from the decay of B^ 

can be approximated by the same expression with the coefficients C^+ '* replaced by C^o '*• 
According to the isospin symmetry relations in Eqs. (75) and (83), this is equivalent to the 

substitutions C^+ ' — )■ C^+ ' 

J/ip TT'^ir" is simply 



and C^+ — )> C^+ . 



Thus the ratio Rn^ for the final state 



1-0+ 



i?o+[J/^7r+7r- 



2(7i - K,ii)i^C 



K+A 



271 z/C 



K+fl 
B+ 



g'^'jiKuC 



B^ 



2(71 






271 z/C 



K+,1 
B+ 



g^lii^uC 



K+,2 
B+ 



(100) 
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This ratio is different from 1 only if C^+ ' ^ C^+ ' , wliicli means tliat tlie 1 = 1 component 
of tlie amplitude for producing the charm mesons at short distance is nonzero. We will 
constrain this ratio to have the value in Eq. (98a). 

If the line shape for [D^D^ti^^D^D^'-)) was completely dominated by the X(3872) reso- 
nance below the D*^D^ threshold, the ratio Rqj^[D^D^t^^ ^D^D^^] would also be given by 
the expression on the right side of Eq. (100). However the measured value of this ratio 
in Eq. (98b) is significantly larger than that for J/-?/; tt'^tt" in Eq. (98a). This is easy to 
understand. Both channels have a peak just below the D*^D^ threshold from the X(3872) 
resonance. However the {D'^D^ir^, D^D^'j) channel also has a threshold enhancement above 
the D*^D^ threshold. At the peak of the threshold enhancement, which is very close to 
the threshold, the ratio in Eq. (100) may be a good approximation. However, above the 
peak in the threshold enhancement, the energy distribution decreases relatively slowly as 
E increases. Because of this high-energy shoulder, there may be a substantial contribution 
from energies where the ratio i?o+ differs significantly from that in Eq. (100). Some idea 
of the energy range over which the threshold enhancement extends can be obtained from 
measurements of the X(3872) resonance in the D*^D^ and D^D*^ decay channels by the 
Belle and Babar Collaborations [15, 16]. The average of their measurements of the po- 
sition is 2.1 ± 1.2 MeV above the D*^D^ threshold. The average of their measurements 
of the width is 3.5]*^}o MeV. These measurements imply that most of the contribution to 
the threshold enhancement comes from the energy region below the D*^D~ threshold at 
z/ii = 8.1 MeV. These measurements of the position and the width in the D*^D^ channel 
should not be interpreted as the actual position and the width of the X(3872) resonance 
for two reasons. First, the X(3872) resonance and the D*^D^ threshold enhancement are 
smeared into a single peak by the experimental resolution. This explains why the measured 
width in the D*^D^ channel is larger than the upper bound in Eq. (54) from measurements 
in the J/ip tt'^tt" channel. This effect also biases the measurement of the position towards 
larger values. Second, in the analysis procedures in Refs. [15, 16], a D^tt^ whose invariant 
mass is close to the mass of the D*^ is constrained to have an invariant mass exactly equal 
to M^Q. This shifts D^D^n^ events just below the D*^D^ threshold to above the threshold, 
which further biases the measurement of the position towards larger values. These two bi- 
ases explain why the measured resonance position in the D*^D^ channel is larger than that 
in Eq. (1) by more than two standard deviations. 

Our prescription for the branching fraction Br+ for the final state {D^D^7T^,D^D^'~f) is 
an integral over E of the appropriate energy distribution: 

Br+[Z^°5°7r°, D^D'-f] = -i- /^''" dE ^[5+ -^ K+ + {D^D%\ D^D^^)], (101) 

r[5+J J_^j,j, dE 

where T[B+] = 4.02 x 10"^° MeV is the total width of the B+. The integrand is the 
contribution to the inclusive energy distribution in Eq. (71) from terms with the factor 
lm.K,{E). We have chosen the integration region somewhat arbitrarily to extend from the 
D^D^-K^ threshold at -5oo = -7.1 MeV to the D*+D- threshold at vu = +8.1 MeV. We 
will constrain the branching fraction in Eq. (101) to have the value in Eq. (99). 

Our prescription for the ratio -Ro+ for {D^D^tt^, D^D^'-f) is the ratio of two integrals like 
the one in Eq. (101): 

n-n n n-n J^;'' dE dr[B' ^ K' + {D^^TT^ , D^D'-f)]/dE 



-5oo 



/_7" dE dT[B+ ^ K+ + (D^DH^, DOD^-f)]/dE 
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The integrand in the numerator differs from that in the denominator by the substitutions 

C^B+ -^ ^B+ an( 
value in Eq. (98b). 



Cg+ ' — )■ C^+ ' and C^+ ' — )■ C^+ ' . We will constrain the ratio in Eq. (102) to have the 



B. B meson decays into K + x'd 

The X(3872) could arise from a fine tuning of the P-wave charmonium state x'd = Xci(2P) 
to near the D*^D^ threshold. However the x'd could also be a separate resonance above 
the D*^D~ threshold. The Zero-Range+Resonance model of Section III D allows for both 
possibilities. If the mass of the x'd i^ sufficiently far above the D*~^D~ threshold, the 
primary effect of the coupling of the x'd ^'^ ^^^ charm mesons is through their contribution 
to its width. Thus information about the decay B -^ K + x'd ^^ ^^^ absence of resonant 
interactions with charm mesons can be used to constrain the coefficients for the B -^ K 
transition. 

In the decay B ^ K + x'd^ ^^e momentum transferred to the mesons in the final state 
is about 1100 MeV. This might be large enough that the decay can be treated using fac- 
torization methods that separate the hard momentum scales comparable to or larger than 
the momentum transfer from the soft hadronic momentum scales. There have been several 
studies of factorization for the decays B ^ K + Xcj, where Xcj = Xcj(l-P), J = 0,1,2, 
is the lowest multiplet of P-wave charmonium states. One might expect the amplitude for 
this process to satisfy a simple factorization formula analogous to that for the decay of a -B 
meson into two light hadrons [45]. In that factorization formula, the factor associated with 
Xcj is the light-front distribution amplitude for the charmonium state. At leading order in 
the relative velocity of the charm quarks, this factor is proportional to /?' (0), the derivative 
of the radial wavefunction at the origin. However, Chao, Song, and collaborators discov- 
ered infrared divergences at next-to-leading order in ag that are proportional to {mc/'mf,y 
[46, 47]. These infrared divergences signal that the factorization formula holds, at best, 
only up to corrections of order (rric/TTif,)^. Beneke and Vernazza studied the factorization of 
the decay amplitude in the limit mb,mc — t- oo with mc/mb fixed [48]. They demonstrated 
that all infrared divergences at next-to-leading order in a^ can be absorbed into the matrix 
element of a color-octet operator. In the asymptotic limit in which Xcj is a Coulombic char- 
monium state, the color-octet matrix element is proportional to R' {0), but in general it is 
an independent nonperturbative factor. Bodwin, Tormo and Lee studied the factorization of 
the decay amplitude in the limit rrif, — ;■ oo with rric fixed [49] . They proved the factorization 
to all orders in ag up to corrections suppressed by mc/mb. Their factorization formula has 
the same structure as that for decays into two light mesons. The factors associated with the 
Xcj are matrix elements of local operators in nonrelativistic QCD. The matrix element at 
leading order in the relative velocity of the charm quark is proportional to i?' (0). 

There have been attempts to calculate the branching fraction for B^ — )■ K^ + x'd- Meng, 
Gao, and Chao used a naive factorization formula with an infrared cutoff on the divergences 
at order ag and obtained the branching fraction 1.8 x 10~^ [50]. They also applied their 
method to the decay into K^ + Xci and obtained essentially the same branching fraction, 
which is smaller by about a factor of 3 than the measured value. Liu and Wang calculated 
the branching fraction for B^ — )■ K^ + x'd [^1] using a perturbative QCD method that 
incorporates Sudakov effects [52]. Their result was (7.9l3'g) x 10~^. They did not calculate 
the corresponding result for K'^ + Xd) so the accuracy of their method can not be judged 
by comparing with data. 
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In the absence of any reliable calculations of the branching fraction for B~^ — ?> K~^ + x'ci^ 
we will use a simple phenomenological estimate. We scale the measured branching fraction 
for B~^ — )■ K~^ + Xci by factors that take into account its dependence on the mass and 
the wavefunction of the charmonium state. The dependence on the mass M^^^ is primarily 
through a multiplicative factor \^^'^(Mb, rnx-, ^xci) ^^^^ comes from the phase-space integral 
and from the Lorentz structure of the decay amplitude. We assume that the dependence on 
the wave function comes primarily from a multiplicative factor |i?' (0)p, where i?' (0) is 
the derivative of the radial wavefunction at the origin. The factor |i?' (0)p for both the 2P 
and IP multiplets has been calculated for four potential models in Ref. [53]. The ratio of 
this factor for 2P and IP ranges from 0.97 to 1.42. We interpret this as a ratio 1.20 with a 
theoretical error ±0.23 that can be added in quadrature with the experimental error in the 
branching fraction for B^ — )■ K^ + Xci- The measured branching fraction for B^ — > K^ + Xci 
is (4.9 ± 0.5) X 10~"^. Multiplying by our two scaling factors, our estimate for the branching 
fraction is 

B,.|B- - A- + x41 - [(5.9 ± L3) x 10"'] ^J^^^^^^^.^^;^ - (103) 

The last factor ranges from 0.53 if the x'ci is near the D*^D^ threshold to 0.39 for M^' = 
4000 MeV. 

Our estimate for the branching fraction in Eq. (103) is based only on information about 
charmonium states. Thus the appropriate limit of the Zero-Range+Resonance model in 
Sec. Ill D is the Flatte limit in which the charm mesons scatter only through their coupling 
to the resonance. In this limit, the partial width for 5+ — )■ K~^ + x'd is given by Eq. (95), 
with the short-distance factor T^ given by Eq. (72). The factor A^/^(Mb, rriK, M^o + ^o) in 
r^ came from taking the limit £" — )■ in \^^'^{Mb, mx, M*o + Mq + E). This factor matches 
the factor \^^'^{MB,mK, M^r ) in Eq. (103). Matching the other factors in Eqs. (95) and 
(103), we obtain the constraint 






*o i=)0 



^.1±1.8) X lO^^l (104) 



C. B meson decays into K + D*"D 



The invariant mass distributions of charm meson pairs in the decays B ^ K + 
[D*^D^, D^D*^) have been measured by the Babar and Belle Collaborations from the D*^D^ 
threshold near 3872 MeV up to 4000 MeV [15, 16]. Measurements of the enhancement 
near the D*^D^ threshold, which is associated with the X(3872), provide the constraints 
on Br+[L)0Z)°7r°,L'°r)°7] and i?o+[^°5°7r°, D^D^] described in Section VIA. The distri- 
butions from 3880 MeV to 4000 MeV provide additional constraints. In particular, they 
constrain the possibility of an additional resonance that could be associated with the x'd- 

The Babar Collaboration used a data sample that consists of about N^b = 3.83 x 10^ 
BB pairs [15]. They combined both B^ and B^ decays and the D^D^-k^ and D^D^j decay 
channels into single measurements in each energy bin. Their efficiencies for B~^ and B'^ 
decays were e+ = 4.3 x 10~^ and cq = 0.7 x 10"*^, respectively. The numbers of events 
were given in energy bins of 2.5 MeV. The Belle Collaboration used a data sample that 
consists of about Nbb = 6.57 x 10^ BB pairs [16]. They combined B~^ and B^ decays 
into separate measurements for the D^D^ti^ and D'^D^'y decay channels in each energy bin. 
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FIG. 1: D*^D^ invariant mass distributions from the decays B -^ K + {D*^D°,D^D*^). The data 
are the numbers of {D*^ D^ , D^ D*^) events per 10 MeV bin observed by the Babar Collaboration 
[15] (open dots) and by the Belle Collaboration [16] (solid squares). The dotted curves are the 
best fits to the data in the bins from 3900 MeV to 4000 MeV with no signal and a background 
proportional to the phase-space volume. 



From the efficiencies that were given, we infer that the efficiencies for B^ and B^ decays 
were e+ = 3.93 x 10"^ and eo = 0.63 x 10"^ in the D°D°7 channel and e+ = 4.24 x 10"^ 
and Co = 0.69 x 10~^ in the D^D'^tt^ channel. The numbers of events were given in energy 
bins of 2 MeV. In the Belle analysis, the energy resolution increased from about 0.5 MeV 
at 3880 MeV to about 1 MeV at 3900 MeV and then to about 2 MeV at 4000 MeV. 

The Babar and Belle measurements provide no evidence for a second resonance above 
3880 MeV, but they also do not exclude it. To simplify our analysis, we combine adjacent 
bins into bins with 10 MeV width. The minimum number of events in any bin is then 
8. This is large enough that the experimental uncertainty in a bin with N events can be 
approximated by viV. The larger bins also decrease the sensitivity to the experimental 
energy resolution. We add the number of Belle events in the D^D^tt^ and D^D^'-f channels 
to obtain the total number of {D'^D^tt^, D^D^'-j) events in each energy bin. The resulting 
Babar and Belle data sets are shown in Figure 1. The curves are the best fits to the data 
from 3900 MeV to 4000 MeV with no signal and with a background proportional to the 
phase-space volume. The large excesses in the ffist bin are associated with the X(3872). 

The predicted number of D*^D^ and D^D*^ events in an energy bin is the sum of the 
predicted numbers of signal events A^sig and background events iVbg. We consider only 
the energy bins above 3880 MeV, so there is a negligible overlap with the X(3872) signal 
region, which is taken into account through the constraints on Br+[D°Z)°7r°, D'^Z)°7] and 
Ro+[D'^D^n^, D^D^j] described in Section VIA. The predicted number of {D^D^tt^, D^D^^i) 
signal events from both B^ decays and B^ decays in an energy bin extending from Ei to E2 
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+eo ^[5° ^ ir° + (Z)°DV0, D^^S)]") . (105) 

We take the background to be incoherent and proportional to the phase-space volumes for 
the two successive 2-body decays B ^ K + {D*D) and {D*D) — > D*^ + -D°. The predicted 
number of background events in an energy bin is 

iVbg = Cbg / dE \'^\Mb, Mk, M,o + Mo + E) \'^\M,o + Mo + E, M,o, Mq), (106) 

J El 

where Cbg is an adjustable constant that is different for the Babar data and for the Belle 
data. 

The inclusive energy distribution dV /dE for B^ decay is given by the factorization formula 
in Eq. (71). In the Zero-Range+Resonance model, the sums are over i,j G {0,1,2}. The 
corresponding energy distribution for B^ decay can be obtained by the interchange C^+ ' 4-)- 

C^+ ' . The energy- dependent functions \v[vFij{E) are expressed as linear combinations of 
the imaginary parts of k{E)^ /^i(-E), 70, 7i, and z/ in Eqs. (42), (84), and (85). The energy 
distribution for [D^D^ti^, D^D^'-y) is the sum of the terms with the factor \m.K{E). The Zero- 
Range-|- Resonance model should give an accurate description of this energy distribution in 
the D*D threshold region. We will use this model all the way up to 4000 MeV, which is 
120 MeV above the D*^D~ threshold. At energies well above the D*^D~ threshold, the 
resonance term in the amplitude may still be accurate, but the zero-range approximation to 
the non-resonance term becomes inadequate. Thus in the high energy region, our model for 
the non-resonance amplitude can at best be regarded as illustrative. 



VII. NUMERICAL ANALYSIS 

In this section, we analyze the available data to determine whether they are able to 
discriminate between binding mechanisms for the X(3872). Our strategy is to fix the in- 
teraction parameters of the Zero-Range-|-Resonance model and then to vary the B ^ K 
transition coefficients to obtain the best possible fit to all the constraints. We repeat this 
procedure for various values of the interaction parameters to see if there are any parameters 
that are preferred by the constraints. 

Given a set of interaction parameters, we determine the B -^ K transition coefficients 
by minimizing the x^ associated with the following constraints: 

• the experimental value for RQ+lJ/ipn^n^] in Eq. (98a), which is set equal to the 
theoretical result in Eq. (100). 

• the experimental value for Ro^[D^D^7T^,D^D^'~f] in Eq. (98b), which is set equal to 
the theoretical result in Eq. (102). 

• the experimental value for Br+fD^D^vr", D^D^'-f] in Eq. (99), which is set equal to the 
theoretical result in Eq. (101). 
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FIG. 2: The niinimuni of x^ ^^ ^ function of i^ for 7 = 28 MeV, g = 0.4, and various values of 71. 

• the constraint in Eq. (104), which follows from our phenomenological estimate for the 
branching fraction for B^ — )■ K^ + x'd- 

• the numbers of B ^ K + [D*^ D^ , D^ D*^) events observed by the Babar and Belle 
collaborations in each of the 12 energy bins covering the range 3880 to 4000 MeV. The 
two data sets are shown in Figure 1. The predicted number of events in an energy bin 
is the sum of the signal given by Eq. (105) and the background given by Eq. (106). 

The total number of constraints is 28. The complex constants C^+ '* are determined up to 
an overall phase by 5 real parameters. The numbers of background events in the Babar and 
Belle data are determined by 2 additional parameters. These 7 parameters are then varied 
to minimize the x^ associated with the 28 constraints. 

In the Zero-Range+Resonance model, the independent interaction parameters are 7, 71, g, 
and u, with 70 given by Eq. (45). The allowed region for the real and imaginary parts of 7 are 
given in Eqs. (57) and (58). Our constraints on the energy distributions are insensitive to the 
value of 7, because they either involve integrals over the D*D threshold region or they involve 
only higher energies. We therefore set the real part of 7 to its central value in Eq. (57) and 
its imaginary part to zero: 7 = 28 MeV. The resonance coupling constant was determined in 
Eq. (59) tohe g = 0.40. The undetermined interaction parameters are the isovector inverse 
scattering length 71 and the resonance parameter z/. Our constraints are insensitive to the 
imaginary parts of u and 71, so we set them to 0. We expect the isovector inverse scattering 
length 71 to be large compared to the momentum scale nu = 125 MeV associated with 
isospin splittings. We will therefore consider values in the region |7i| > 2k,ii = 250 MeV. 
We consider values of i' in the range —10 MeV < Keu < 140 MeV. For u < —10 MeV, there 
would be a second resonance below the D*^D^ threshold in addition to the X(3872), which 
is not observed. For u > 140 MeV, the second resonance would be above 4000 MeV, beyond 
the region covered by the Belle and Babar data. 

We have performed fits of the 7 parameters to the 28 constraints for 7 = 28 MeV, g = 0.4, 
and different values of z/ and 71. The parameter u has been varied from —10 MeV up to 140 
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MeV in steps of 1 MeV. Since our results are less sensitive to 71, we have repeated the fit 
only for a few discrete values of 71: ±250, ±500 and ±1000 MeV. The minimum values of x^ 
for the various values of 71 are shown in Figure 2 as a function of v. For v < 100 MeV, the 
minimum x^ lies in the range 30 < Xmm < 39 for all values of 71, indicating that the quality 
of the fit is rather insensitive to 71 and u. The minimum value of Xmin ~ 21.7 is reached at 
z/ = 119 MeV and 71 = —1000 MeV. For these values of z/ and 71, we verified that the fit 
and the value of Xmm ^^^ insensitive to variations in the real and imaginary parts of 7. For 
u > 130 MeV, the minimum x^ lies in the range 29 < Xmm < 36, indicating again that the 
quality of the fit is rather insensitive to 71. 

In Figure 3, we show the D*^D^ invariant mass distributions for 7 = 28 MeV, g = 0.4, 
and three pairs of parameters (z^, 7i) that are approximate local minima of x^'- (~10,250), 
(50, —250), and (120, —1000) MeV. The third pair of parameters is close to the global mini- 
mum of x^- The predicted numbers of D*^D^ events in a 10 MeV bin observed by the Belle 
(Babar) experiment can be obtained by averaging the upper (lower) solid line over the bin. 
The background contributions are shown as dotted lines. For the first set of parameters, the 
X(3872) arises from a fine-tuning of the energy of the x^^. For the second and third sets of 
parameters, the X(3872) arises from a fine-tuning of the strength of the interactions between 
D*^ and /)". For all three pairs of parameters, there is an X(3872) resonance below the 
D*^D^ threshold and a D*^D^ threshold enhancement above the threshold. For the first set 
of parameters, there are no additional peaks in the energy distribution. For the second set of 
parameters, there is an additional peak near 3913 MeV corresponding to the x'd resonance. 
For the third set of parameters, there is a peak near 3983 MeV corresponding to the Xci 
resonance, but it is preceded by a sharp decrease in the number of events down to near the 
background level. This surprising behavior arises from destructive interference between the 
resonant and nonresonant amplitudes for producing charm meson pairs. 

In Figure 4, we show the predicted energy distributions for the {D^D^n'^, D^D^'j) 
signal and for the {D^D^n^ , D^D^n^, D^D^n^, D^D~'j) signal from both B^ de- 
cay and B^ decay for the same three sets of parameters as in Figure 3. The 
{D'^D^n" ,D'^D~7i~^,D^D~'K^,D'^D~'~f) signal is the sum of the terms in the factorization 
formula with a factor of lm.K,i{E). The two signals correspond to {D*^ D^ , D^ D*^) and to 
(D*^ D^ , D^ D*^) above the appropriate thresholds. For all three sets of parameters, the 
energy distribution for {D^D^n'^ , D^D~7t~^ , D^D^n^, D^D^'-f) has an X(3872) resonance 
below the D*^D^ threshold, but there is no D*^D^ threshold enhancement just above the 
threshold. The resonance should be interpreted as a contribution from D^D~'y, because it 
is below the thresholds for D~^D^tt^ and D~^D~7r^. A quantitative treatment of this con- 
tribution would require using the energy-dependent width for the D*~^ defined in Ref. [9]. 
For the second set of parameters, there is a peak at the x'ci resonance near 3913 MeV in 
both D*^D^ and D*~^D~ , but the peak is higher in D*^D~ . For the third set of parameters, 
there is constructive interference in D*^D~ at the x^j^ resonance near 3980 MeV in contrast 
to the destructive interference in D*^D^. 

In Figure 5, we show the behavior of the line shapes in Figure 4 in the D*D threshold 
region. At the X(3872) resonance, the line shapes in [D^D^ii^ ^D^D^'-y) and in D~^D~'y 
coming from B~^ decay and 5° decay differ only in their normalizations. The line shapes in 
{D'^D^TT^, D^D^'j) have a D*^D^ threshold enhancement whose peak is about as far above 
the D*^D^ threshold as the X(3872) resonance is below the threshold. The two maxima 
are well-separated, because we have set Im7 = 0. If Im7 is large enough, the resonance and 
the threshold enhancement will be smeared into a line shape with a single maximum. The 
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FIG. 3: Numbers of B ^ K + (D*^ D^ , D^ D*^) events as functions of the invariant mass of the 
charm meson pair. The data are the numbers of events observed by the Babar Cohaboration [15] 
(open circles) and by the Belle Collaboration [16] (solid squares) in 10 MeV bins. The predicted 
{D^D^7r^,D^D^j) distributions for Belle (upper solid lines) and for Babar (lower solid lines) are 
shown for 7 = 28 MeV, g = 0.4, and three sets of values of (i/, 71). The background contributions 
are shown as dotted lines. 
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FIG. 4: Energy distributions for the {D^D^n^, D^D^'y) signal (solid lines) and for the 
{D~^D^7T~ ,D^D^7:^, D^D^7r^,D~^D~'y) signal (dashed lines) from both i?+ decay (thicker lines) 
and B^ decays (thinner lines) for 7 = 28 MeV, g = 0.4, and three sets of values of (z/, 71). 
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FIG. 5: Line shapes in the D*D threshold region for the {D^D^7r^,D^D^^) signal (solid lines) and 
for the {D^ D^7r~ , D^ D~ TT^ , D~^ D~ TT^ , D^D^^) signal (dashed lines) from both B~^ decay (thicker 
lines) and B^ decays (thinner lines) for 7 = 28 MeV, g = 0.4, and three sets of values of (2^,71). 
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lines shapes from B~^ decay and B^ decay can be quite different in the region up to D*~^D 
threshold, and they can change dramatically at that threshold. The energy distributions in 
{D+D^n-,D^D-n+,D+D-7r^,D+D--f) are close to zero between the D*^D° and D*+D- 
thresholds, but they can increase dramatically above the D*^D^ threshold. 

VIII. SUMMARY 

We have derived general line shapes for the X(3872) that can be used to discriminate 
between alternative binding mechanisms. In the D*^D^ threshold region, which extends 
only to a few MeV from the threshold, the line shapes are universal in the sense that they 
are determined only by the binding energy and width of the X(3872) [10]. In the D*D 
threshold region, which extends to tens of MeV from the threshold, the line shapes depend 
on the binding mechanism. The mechanism could be the fine-tuning of the interaction 
between the charm mesons to near the critical strength for a bound state, in which case the 
appropriate line shapes are the zero-range line shapes of Ref. [10]. The mechanism could also 
be the fine-tuning of the energy of a resonance to near the D*^D^ threshold, in which case 
the appropriate line shapes are the Flatte line shapes. The Flatte line shapes are essentially 
those introduced in Ref. [20], except that one cannot ignore the essential contribution to 
the line shape in the {D^D^tt^, D^D^'-f) channel from the X(3872) resonance peak below the 
threshold. Our general line shapes include the zero-range line shapes and the Flatte line 
shapes as special cases. 

Our general line shapes take into account scattering between two coupled channels with 
the same reduced mass /i. This is a good approximation if the masses of the pairs of particles 
in the two channels differ only by isospin splittings. These line shapes could be applied to any 
weakly-bound hadronic molecule that has S-wave couplings to a pair of thresholds separated 
by isospin splittings. However the complications associated with the coupled channels are 
only relevant if the binding energy of the molecule and the widths of the constituents are 
all much smaller than the isospin splittings, which are typically less than 10 MeV. The 
X(3872) may be the unique hadron that satisfies this requirements. If these conditions 
are not satisfied, one might as well ignore the complications associated with the isospin 
splittings. The zero-range approximations for the direct interaction between the pair of 
mesons may also be inadequate in this case. 

A crucial ingredient in our general line shapes are the coupled-channel scattering ampli- 
tudes for the neutral and charged charm meson channels [D*^D'^)+ and (D*+D^)+, which 
are given in Eqs. (39). The resonance propagator is given in Eq. (41). They depend on 4 
interaction parameters: 70, 71, z/, and g. The scattering amplitudes satisfy the constraints of 
unitarity exactly if these parameters are real. The analytic continuation of the parameters 
to complex values can be used to take into account some of the effects of states that are not 
treated explicitly. The line shapes in the universal region depend only on the inverse scatter- 
ing length 7, which satisfies Eq. (43). It is therefore convenient to use Eq. (45) to eliminate 
7o in favor of 7. The real and imaginary parts of 7 are constrained by measurements of the 
binding energy and width of the A(3872). The constraints are given in Eqs. (57) and (58). 
We expect the isovector inverse scattering length 71 to be much larger than /ti ^ 125 MeV, 
but we do not have any useful quantitative constraints on this parameter. If the resonance 
is identified with the P-wave charmonium state x'ci — Xci(2-P), charmonium phenomenol- 
ogy can be used to constrain the resonance parameters u and g. The coupling constant g 
is determined from phenomenological models of the decays x'ci ~^ D*D to have the value 
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g = 0.4 given in Eq. (59). The variations among potential models are sufficiently large that 
they do not provide any useful constraints on the energy parameter u. 

The general line shapes for a weakly-bound hadronic molecule that is produced by a 
short-distance process can be expressed in terms of a factorization formula for the inclusive 
energy distribution. In the case of the X(3872), an example of a short-distance process 
is a B ^ K transition with momentum transfer of about 1100 MeV, which can produce 
pairs of charm mesons in the D*D threshold region. The factorization formula for the 
inclusive energy distribution for states with quantum numbers l"*""*" that are produced by 
the B~^ — )■ K~^ transition is given in Eq. (71). All the dependence on the energy E is in the 

long-distance factors lm.Fij{E). The short-distance coefficients C^+ '* are complex constants. 
The factorization formula for another short-distance process differs only in the values of the 
short- distance coefficients. The short-distance coefficients for the B^ — )■ K'^ transition are 
determined in terms of those for the B~^ — )■ K~^ transition by the isospin-symmetry relations 
in Eqs. (75) and (83). These short-distance coefficients are constrained by measurements 
oi B -^ K + X(3872), such as those in Eqs. (98) and (99). The additional constraint in 
Eq. (104) was obtained from a phenomenological estimate of the decay rate for B — )■ K + x'ci- 

In the factorization formula for the inclusive energy distribution in Eq. (71), the long- 
distance functions ImFij {E) can be resolved into contributions from different decay channels 
of the X(3872) by expressing them in forms consistent with the Cutkosky cutting rules. The 
functions lm.Fij{E) are expressed as linear combinations of the imaginary parts of the func- 
tions n{E) and ni{E) and the imaginary parts of the parameters 70, 71, and v in Eqs. (42), 
(84), and (85). The terms proportional to \m.K{E) can be interpreted as the contributions 
from {D^D^7T^,D'^D^'~f). The terms proportional to lm.K,i{E) can be interpreted as the con- 
tributions from {D'^D^TT' ,D^D^ 71^,0^ D^ 71^, D~^D~'~f). The terms proportional to Imz/, 
Im7o, and Im7i can be interpreted as the contributions from other decay modes of x'cn 
other channels with isospin 0, and other channels with isospin 1. 

In Section VII, we used our line shapes to carry out a phenomenological analysis of the 
data on B ^ K + X(3872) and data from the Belle and Babar Collaborations on 5 — )■ 
K + D*^D^, with the D*^D^ invariant mass ranging up to 4000 MeV. We assumed that the 
binding mechanism for the X(3872) is either a fine-tuning of the strength of the interaction 
between the charm mesons or a fine-tuning of the energy of the P-wave charmonium state 
Xci- We extrapolated our line shapes for the D*D threshold region all the way up to 128 MeV 
above the D*^D^ threshold. In the high energy region, the nonresonant contributions to the 
production amplitudes, which were derived using a zero-range approximation, can at best 
be regarded as an illustrative model. We found that the two mechanisms for the binding 
of the X(3872) are both compatible with the data for B -^ K + D*^D^ and our other 
constraints. One way to exclude the tuning of the Xci energy as a binding mechanism 
for the X(3872) is to observe the Xci resonance as a separate peak in the invariant mass 
distribution for D*^D^ and D*^D~ . In the Belle and Babar data, there is no obvious peak in 
the D*^D^ invariant mass distribution between 3880 MeV and 4000 MeV. However, the width 
of this peak is determined by the resonance parameters z/ and g and its height is determined 
by the short-distance coefficients C^+ '*. Our calculations show that a x'ci resonance in 
this region is compatible with the data for B ^ K + D*^D^ and with our constraints on 
the short-distance coefficients. An interesting possibility for this x'd resonance is to have 
destructive interference between the resonant and nonresonant amplitudes in the D*^D^ 
decay channel and constructive interference in the D*^D^ decay channel. This possibility 
is actually realized for the parameters that minimize the x^ for our constraints, which give 
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a x'ci resonance near 3980 MeV. 

An alternative way to exclude the x'd resonance mechanism for the binding of the X(3872) 
is to calculate the x'd mass using lattice gauge theory. The masses of excited charmonium 
states and states in the cc meson spectrum with exotic quantum numbers have been calcu- 
lated by Dudek et al. using lattice gauge theory without dynamical light quarks [54]. The 
masses of the 2P charmonium multiplet are significantly higher than those of the ID multi- 
plet. Their analysis suggests that the 2"'""'" cc meson discovered near 3930 MeV is more likely 
to be the ^D2 ground state than the first radial excitation of the ^P2- This suggests that the 
2P multiplet, including x'ci^ has higher mass. If these results are confirmed by lattice calcu- 
lations with dynamical light quarks, it would exclude the x'd resonance mechanism. Other 
resonance mechanisms, such as the tuning of the energy of a 1"'""'" tetraquark cc meson, are 
not easily constrained by phenomenology, but they can also be ultimately ruled out using 
lattice QCD calculations. 

We used charmonium phenomenology to constrain the resonance parameters g and u. 
Phenomenological estimates for the scattering parameters 70 and 71 could be obtained from 
meson potential models, such as those in Refs. [25, 38-42], or from meson scattering models, 
such as the one in Ref. [43]. Thus far, these models have been used primarily to calculate 
binding energies. They could also be used to calculate scattering variables, such as 70 and 
7i. A particularly convenient pair of scattering variables are the residues Zq and Zi of the 
poles in the elastic scattering amplitudes foo{E) for D*^D^ and fu{E) for D*~^D~ at the 

X(3872) resonance. For the Zero-Range+Resonance model, the ratio Z-^^ /Z^ is given in 
Eq. (47). It does not depend on the resonance parameters g and u, so it is determined 
primarily by 71. 

Our general line shapes could be used by experimentalists to carry out a global analysis 
of the energy distributions for the decay modes of the X(3872) that is not biased towards 
a specific binding mechanism. One complication is the number of independent parameters 
in the line shapes. The line shapes in [D^D'^tt'^, D^D^'-j) depend essentially on 5 interaction 
parameters: the real parameters 71, g, and v and the complex parameter 7. For other decay 
channels, there is also a normalization parameter, such as (Im7i)'^/'^'^ '^ for J/iIj-k^-k" and 
(Im7o)"'/'^'^ '^ ^ for J / tp n^ ti~ n'^ . The imaginary parts of 70, 71, and u give contributions 
to Im7, as indicated by Eq. (46), but their effects should otherwise be negligible in the D*D 
threshold region. In addition to the interaction parameters, there are 5 real parameters 
associated with the short- distance coefficients for the production of the charm meson pairs 
and the resonance. In an analysis of the line shapes of the X(3872) produced by B decays, 
the determination of the 5+ interaction parameters and the 5 short-distance parameters 
would require analyzing several decay channels in both B^ and B^ decays. 

An alternative strategy would be to carry out two separate analyses assuming either the 
dynamical mechanism or the resonance mechanism for the binding of the X(3872). With the 
dynamical mechanism, one could use the zero-range line shapes, which have 3+ interaction 
parameters and 3 short- distance parameters. With the resonance mechanism, one could use 
the Flatte line shapes, which have 3+ interaction parameters and 1 short-distance parameter. 
If one set of line shapes gives a significantly better global fit to the data, it would be evidence 
in favor of the corresponding binding mechanism for the X(3872). 
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Appendix: Renormalization of the Zero-Range+Resonance Model 

The Zero-Range+Resonance model considered in Sections III D and V E can be derived 
from a renormalizable quantum field theory with local interactions. In this Appendix, we 
consider the renormalization of the interactions of this field theory and the renormalization 
of the local operators that describe production by a short-distance process. 

A. Renormalization of the Interactions 

The Zero-Range+Resonance model can be represented by a nonrelativistic field theory. 
The fields are scalar fields for the spin-0 charm mesons D'^, D^, D~^, and D^ , vector fields 
for the spin-1 charm mesons D*^, Z)*°, D*~^, and D*~ , and a vector field for the resonance 
X- The vector fields, such as x™"? have a Cartesian vector index m. The interaction terms 
in the Hamiltonian density are 



+ ^0 X'"^X". 



271 - - I2-JT 

H,^, = —(D*Dr^Ao(D*Dr + \ — \(D*Dr^Gox" + x"^ G^ {D* 0)""] 
/i V /^ 

(107) 



where [D*D)"^ is a 2-component column vector with a Cartesian vector index m. Its upper 
and lower entries are the combinations of fields that annihilate pairs of charm mesons in the 
neutral channel {D*'^D^)^ and in the charged channel {D*'^D^)^ given in Eqs. (2). The bare 
parameters are the three independent entries of the symmetric matrix Aq, the two entries of 
the column vector Gq, and uq. The factors of y^27r/fi in the Hamiltonian density have been 
inserted to simplify the renormalization equations. The x^^X™ term has been included in 
the interaction Hamiltonian, because its coefficient z/q requires renormalization. 

The matrix of transition amplitudes A{E) for charm meson pairs in the channels 
(iI)*°Z)°)_|_ and {D*~^D~)^ can be calculated by summing diagrams constructed out of the 
two types of vertices, resonance propagators, and loop subdiagrams involving pairs of charm 
mesons. The ultraviolet divergences in the loop diagrams can be regularized with an ul- 
traviolet momentum cutoff Auv- The matrix A{E) can also be derived by solving the 
Lippmann-Schwinger integral equation, as described in Ref. [55]. It is convenient to express 
the solution in terms of the matrix of scattering amplitudes f{E) defined by Eq. (12), which 
differs from A{E) by a factor of /i/(27r). The solution for the inverse of f{E) is 

fiE)-' = - f Ao + Go-^^G^) - -Auv/ + K{E), (108) 

V E - Vq J TT 

where / is the identity matrix and K{E) is the diagonal matrix given in Eq. (14). The 
energy dependence on the right side of Eq. (108) allows f{E)^^ to be expressed in the form in 
Eq. (35), in which the dependence on Auv has been absorbed into renormalized parameters. 
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The renormalized parameters are the three independent entries of the symmetric matrix A, 
the two entries of the column vector G, and v. The relations between the renormalized 
parameters and the bare parameters are given by 



Z-^Ao, 
Z Go, 

where the renormalization matrix Z is 



A 
G 

V 



Z-') Ao'Go, 



Z = 1+ -AuvAq. 



vr 



(109a) 
(109b) 
(109c) 



;iio) 



The verification of the equality between the expressions for f{E)~^ in Eqs. (108) and (35) can 
be simplified by taking advantage of the existence of renormalization-invariant combinations 
of parameters. Rewriting Eq. (109a) as Z~^ = AA(^^ and inserting it into Eqs. (109b) and 
(109c), we find that the following combinations of parameters are renormalization invariants: 



Ao'Go 
GqAq 

Gr, An G(] 



A'^G, 

G^A-\ 

V - G^A-^G. 



^0 ~ "-^o^^o "-^0 
The following function of E is also a renormalization invariant: 



An + Gr 



E -vq 



Gq 



Ao^ 



A + G 



E 



-G 



■T 



-1 



A' 



(Ilia) 
(111b) 
(111c) 



;ii2) 



This can be verified by multiplying both sides on the left by Aq + GqGq /{E — z/q) and on 
the right by A + GG'^ /{E — v) and then using the renormalization invariants in Eqs. (111). 
The equality between the expressions for f{E)~^ in Eqs. (108) and (35) follows immediately 
from Eq. (112). 

The complete propagator for the resonance x can be obtained by summing the geometric 
series of self-energy corrections. The self-energy can be obtained by summing a geometric 
series of one-loop diagrams: 



E(^) = -G^ { -Auv/ - K{E) 

vr 



I + Aoi-AvyI-K{E) 



IT 



-, -1 



Gq. 



(113) 



The complete resonance propagator is therefore 

'2 



P{E) 



E -uq 



Gq —AuyI 

vr 



K{E) A 



^Auv/ - K{Ey\ A,'Go 



This can be expressed in the form 

P{E) =\e-uo + G^A,'Go - G^Ao 1 (A-^ - KiE))-' A.'Go 



-1 



;ii4) 



(115) 



By using the renormalization invariants in Eqs. (HI), this propagator can be expressed 
in terms of renormalized parameters and then simplified to the renormalized expression in 
Eq. (36). 
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B. Renormalization of the Production Operators 

The production of particles at short distances can be represented in an effective field 
theory by local operators that create the particles when acting on the vacuum. The ma- 
trix elements of the local operators are in general ultraviolet divergent. The corresponding 
renormalized operators are linear combinations that have finite matrix elements. The renor- 
malized operators can be determined by calculating matrix elements of the local operators 
and then constructing linear combinations whose matrix elements are finite. Alternatively, 
the renormalized operators can be deduced by inspired guesswork. 

In the Zero-Range+Resonance model, the particle that are produced at short distances 
can be pairs of charm mesons or the resonance. There are three leading operators: the two 
components of the column vector (D*/))™ that appears in the interaction Hamiltonian in 
Eq. (107) and the resonance field x"*- They all have a Cartesian vector index m. We denote 
the corresponding renormalized operators by (9™, O^, and O^. The Green's function Fij{E) 
for pairs of the operators O^ are defined in Eq. (70). The Green's function for a pair of 
operators x™ is just —P{E), where P{E) is the resonance propagator in Eq. (114). Since 
this propagator can be expressed in the renormalized form in Eq. (36), the resonance field 
has finite matrix elements. We can therefore choose its hermitian conjugate as one of the 
renormalized operators: 

Of = X™. (116) 

The matrix elements of [D*D)"^ are ultraviolet divergent. The renormalized operators O^ 
and O™ are linear combinations of the hermitian conjugates of all three local operators. 
A particular convenient choice consists of the operators obtained by differentiating Tiint 
with respect to the components of {D*D)"^. The Green's functions for these operators 
are proportional to the scattering amplitudes. Our choices for the last two renormalized 
operators are 



Of\ [2^ 



'0 

■)mt 



or J V /i 



Ao{D*Dr + Gox'^. (117) 



The normalization factor has been chosen so that the Green's functions for these operators 
are exactly equal to the scattering amplitudes fij{E) for i,j G {0, 1}. The remaining Green's 
function Fi2{E) for i = 0,1 are then given by 

Fi,{E) = [f{E) [{E - z/o)Ao + G.Glr^ G,)^ , (118) 

where f{E) is the matrix of scattering amplitudes. The corresponding renormalized ex- 
pression is obtained by replacing z/q, Aq, and Gq by u, A, and G. The equality of the two 
expressions for Fi2{E) follows from the fact that Aq^Gq = A^^G is an eigenvector of both 
z/qJ - Aq^GoG^ and ul - A'^GG'^: 

(z/o/ - A,^G,Gl) A^'Go = (z/o - G^A.'Go) A.'Go, (119a) 

{ul - A-^GG'^) A-^G = {u-G^A-^G)A-^G. (119b) 

The Green's functions Fij{E) for i,j G {0,1} are not exactly equal to the scattering 
amplitudes fij{E). The diagrams contributing to the Fij{E) are the same as the diagrams 
for fij{E) except that the leading order diagram for scattering of charm meson pairs is 
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omitted. Thus Fij{E) actually differs from fij{E) by the additive constant (Ao)ij. This 
constant gives an energy-independent contribution to the inclusive energy distribution in 
Eq. (71). This additive contribution can be interpreted as a constant background in decay 
channels of the X(3872) that are not treated explicitly. 

The renormalized operators (9™ and O^ defined by Eq. (117) have different dimensions 
from the renormalized operator O^ = x"^- For some purposes, it is more convenient for all 
the renormalized operators to have the same dimensions. One way to arrange for (9™ to have 
the same dimension as (9™ and O^ is to multiply x"^ by a constant with the appropriate 
dimension. In the text, the multiplying factor is chosen to be (G'"^G'/2)^/^, where G is 
the two-component vector of renormalized interaction parameters. With this choice, the 
renormalized expressions for F22{E) and Fi2{E) are given by Eqs. (80) and (81). 
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